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Abstract 

APPLICABILITY  OF  THRUST  AUGMENTING 
EJECTORS  IN  A  SUPERSONIC  CRUISE 
CONFIGURATION 

by  Wendell  Scott  Hertzelle 

Chairperson  of  the  Supervisory  Committee:  Professor  Abraham  Hertzberg 
Department  of  Aeronautics  and  Astronautics 

Thrust  augmenting  ejectors  were  analyzed  by  the  author  in  an  attempt  to  see  if  they  could 
be  used  to  provide  a  significant  thrust  increase  over  the  baseline  thrust  of  a  primary  core. 
This  was  done  for  the  purpose  of  determining  if  leaving  the  ejectors  open  in  the  cruise 
configuration  of  the  High  Speed  Civil  Transport  would  lead  to  a  thrust  increase  that 
would  at  least  offset  the  weight  of  the  ejectors  themselves.  If  this  was  found  to  be  true, 
then  the  fuel  economy  of  the  HSCT  might  be  improved  by  leaving  them  open  and  not 
closing  them  during  cruise. 

In  analyzing  the  ejectors,  no  assumptions  were  made  regarding  inlet  nor  outlet 
configurations,  so  an  attempt  was  made  to  find  the  point  of  optimal  thrust  augmentation 
by  varying  secondary  stream  bypass  Mach  number  and  the  amovmt  of  flow  entrainment. 
Two  solutions  were  found  to  each  mixing  scenario,  one  subsonic  and  the  other 
supersonic.  These  two  solutions  were  each  analyzed  and  ones  not  satisfying  the  Second 
Law  of  Thermodynamics  were  eliminated.  Analytic  diffuser  and  bleed  losses  were  also 
explored  in  the  analysis  of  the  ejector  flow.  Within  the  limitations  of  the  assumptions 
discussed  in  this  paper,  appreciable  thrust  augmentations  have  been  discovered  over  a 
large  range  of  bypass  Mach  numbers  and  entrained  mass  flows.  This  lead  the  author  to 
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the  conclusion  that  ejectors  warrant  further  research  beyond  a  first  order  analysis,  and 
serious  thought  should  be  put  into  leaving  them  open  in  cruise. 
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CONFIGURATION 
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Chairperson  of  the  Supervisory’  Committee:  Professor  Abraham  Hertzberg 
Department  of  Aeronautics  and  Astronautics 

Thrust  augmenting  ejectors  were  analyzed  by  the  author  in  an  attempt  to  see  if  they  could 
be  used  to  provide  a  significant  thrust  increase  over  the  baseline  thrust  of  a  primary'  core. 
This  was  done  for  the  purpose  of  determining  if  leaving  the  ejectors  open  in  the  cruise 
configuration  of  the  High  Speed  Civil  Transport  would  lead  to  a  thrust  increase  that 
would  at  least  offset  the  weight  of  the  ejectors  themselves.  If  this  was  found  to  be  true, 
then  the  fuel  economy  of  the  HSCT  might  be  improved  by  leaving  them  open  and  not 
closing  them  during  cruise. 

In  analyzing  the  ejectors,  no  assumptions  were  made  regzu-ding  inlet  nor  outlet 
configurations,  so  an  attempt  was  made  to  find  the  point  of  optimal  thrust  augmentation 
by  varying  secondary  stream  bypass  Mach  number  and  the  amount  of  flow  entrainment. 
Two  solutions  were  found  to  each  mixing  scenario,  one  subsonic  and  the  other 
supersonic.  These  two  solutions  were  each  analyzed  and  ones  not  satisfying  the  Second 
Law  of  Thermodynamics  were  eliminated.  Analytic  diffuser  and  bleed  losses  were  also 
explored  in  the  analysis  of  the  ejector  flow.  Within  the  limitations  of  the  assumptions 
discussed  in  this  paper,  appreciable  thrust  augmentations  have  been  discovered  over  a 
large  range  of  bypass  Mach  numbers  and  entrained  mass  flows.  This  lead  the  author  to 


the  conclusion  that  ejectors  warrant  further  research  beyond  a  first  order  analysis,  and 
serious  thought  should  be  put  into  leaving  them  open  in  cruise. 
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PREFACE 


During  the  fall  of  1995.  while  taking  an  "self-structured”  course  in  fluid  mechanics  at  the 
University  of  Washington,  the  author  and  his  class  came  upon  the  topic  of  ejectors  and 
how  they  could  be  used  for  thrust  augmentation.  This  topic  was  discussed  in  class  as  a 
comparison  to  ducted  rockets,  ramjets  and  other  propulsion  topics,  but  it  was  never 
applied  directly  to  the  topic  of  use  on  the  High  Speed  Civil  Transport. 

The  course  advisor.  Professor  Abraham  Hertzberg  mentioned  to  the  author  that  he  had 
been  working  with  another  graduate  student,  Mr.  Greg  Williams,  on  the  use  of  ejectors  to 
increase  thrust  on  the  HSCT.  According  to  Boeing,  ejectors  were  being  looked  at  for  use 
in  noise  abatement  during  takeoff.  By  mixing  the  hot  primary  flow  with  a  cooler 
secondary  flow  inside  the  ejector,  the  noise  that  resulted  from  that  mixing  could  then  be 
absorbed  by  noise  abating  material  contained  in  the  ejector  itself,  thus  reducing  takeoff 
noise. 

Professor  Hertzberg  indicated  that  Boeing  planned  to  leave  the  ejectors  open  during 
takeoff  and  the  subsonic  portion  of  the  takeoff  but  then  close  the  ejectors  as  cruise  was 
reached.  However,  he  felt  that  it  would  be  worthwhile  to  consider  leaving  them  open 
during  cruise  and  see  if  some  thrust  augmentation  could  be  realized.  By  doing  so,  some  if 
not  all  of  the  weight  of  the  ejectors  themselves  would  be  offset  and  therefore  improve  the 
fuel  economy  and  the  range  of  the  aircraft  itself. 

Boeing  was  asked  for  data  regarding  the  cruise  configuration,  the  primary  core  engine 
data,  and  the  ejector  itself.  Taking  this  data,  Greg  Williams  compiled  an  initial 
spreadsheet  that  modeled  the  use  of  the  ejectors  in  the  supersonic  cruise  configuration  at 


IX 


Mach  2.4.  and  from  this  point,  the  author  of  this  paper  took  over  and  expanded  upon  his 
research. 

Modifications  were  made  to  the  baseline  spreadsheets  to  include  diffuser  losses,  bleed 
losses  and  to  calculate  entropy  changes.  Also,  further  thought  was  given  to  the 
interpretation  of  the  mixing  results.  All  of  these  additions  culminated  in  the  creation  of 
this  paper  and  the  results  presented  within. 
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INTRODUCTION 


In  order  to  provide  an  improvement  in  the  design  of  the  High  Speed  Civil  Transport 
(HSCT).  ejectors  were  analyzed  to  determine  if  they  could  be  used  for  thrust 
augmentation.  During  takeoff,  ejectors  were  planned  to  be  used  for  purposes  of  noise 
abatement,  and  then  closed  during  cruise.  How'ever,  the  author  and  his  advisor  wanted  to 
look  at  what  would  occur  if  they  were  left  open  during  the  duration  of  the  flight. 

The  choice  of  cruise  location  and  speed  were  set  by  the  design  specifications  provided  by 
the  Boeing  Corporation.  Their  design,  which  was  already  in  progress,  provided  the 
starting  information  needed  to  begin  the  problem  analysis.  After  being  provided  data  by 
Boeing,  the  author  was  able  to  begin  his  analysis  of  the  applicability  of  ejectors  for  thrust 
augmentation. 


CHAPTER  1 :  BACKGROUND  AND  PROBLEM  SETUP 


1.1:  PROBLEM  ORIGIN 

In  order  to  design  a  supersonic  commercial  transport,  the  aircraft  must  conform  to  a  rigid 
noise  regulation  standard.  FAR  Part  36.  Stage  111.  The  Boeing  Company  has  proposed  to 
use  ejectors  on  their  High  Speed  Civil  Transport  (HSCT)  for  purposes  of  noise 
suppression.  The  phenomena  of  noise  suppression  will  not  be  discussed  in  this  thesis; 
however,  the  thrust  augmentation  properties  of  ejectors  will  be  covered. 

Boeing  plans  to  stow  the  ejectors  in  cruise  configuration,  but  through  analysis  it  can  be 
shown  that  a  significant  increase  can  be  achieved  over  the  baseline  thrust  if  the  ejectors 
are  not  stowed.  The  thrust  increase  could  be  used  to  offset  the  weight  of  the  ejector  itself 
(which  Boeing  stated  is  on  the  same  order  as  the  engine's  weight). 

1.2:  SOLVING  THE  PROBLEM 

1.2.1:  Boeing  Data 

In  order  to  analyze  the  thrust  augmentation  properties  of  ejectors,  the  engine  that  Boeing 
proposes  to  use  for  the  HSCT  was  used.'  According  to  the  data  provided  about  the 
engine,  it  would  have  the  following  fi.xed  characteristics: 
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Table  1 :  Engine  Characteristics  as  Supplied  by 
Boeing 


•  • 

Mass  Flow;  m  ~  198.6  kg/s  =  mx 

Stagnation  Temperature,  Tq  =  1265  Y^=To4 
Stagnation  Pressure,  Po  =  323,780  Pa  =  Po4 

Also,  the  analysis  was  conducted  at  50,000  feet  (15,240  m)  in  standard  atmosphere,  with 
a  flight  Mach  number  of  2.4.  Values  for  ambient  conditions  were  given  in  the  Boeing 
data,  but  Introduction  to  Flight  by  Anderson  was  referenced  to  check  the  accuracy  of  their 
values.^  So,  in  summary,  the  free  stream  conditions  are: 


Table  2:  Free  Stream  Conditions  as  Supplied  by 
Boeing 


Altitude  =  50,000  ft  =  15240  m 
Density,  =  0. 1 865  kg/m^  s  pj 
Ambient  Temperature,  Ty,=  2\6.1  ¥^  =  Ti 
Ambient  Pressure,  =  11 ,597  Va.  =  Pi 
Flight  Mach  Number,  M:c=  2.4  s  My 

1.2.2:  Problem  Setup 

In  order  to  analyze  the  primary  flow  and  the  entrained  secondary  mass  flow,  the  following 
stations  were  defined: 

1 .  Ambient  Conditions 

2.  Just  prior  to  Diffuser  Throat 

3.  Just  after  the  Diffuser  Throat 


4.  Primaiy  Stream  Flow 


5.  Ejector  (Secondary  Stream)  Flow 


6.  Mixing  Plane 

7.  Mixed  Flow^  expanded  to  Ambient  Conditions 


©CD  O 


Figure  1:  Diagram  of  Station  Numbers 

These  stations  were  initially  designated  by  Mr.  Greg  Williams,  a  former  graduate  student 
at  the  University  of  Washington,  when  he  created  the  initial  spreadsheet  that  w'as  used  to 
analyze  this  problem.  So,  to  maintain  consistency  with  previous  work,  the  same 
convention  has  been  adopted. 

To  simplify  the  analysis,  the  following  assumptions  were  made: 

1 .  There  is  no  shock  in  the  secondary'  stream  diffuser. 

2.  The  primary  and  secondary'  flows  are  perfectly  mixed  at  station  6. 

3.  The  static  pressures  at  stations  4  and  5  are  equal  (i.e.  the  primary'  flow,  station  4, 
completely  expands). 

4.  Perfect  gas  relations  hold,  and  all  fluids  are  compressible. 

5.  The  specific  heat  and  specific  heat  ratios  for  the  combustion  products  in  stream  4 
are  identical  to  air.  so  y  =  1 .4  and  R  =  287.0  J/(kg  '  K). 
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6.  Air  is  inviscid. 

7.  Flow  is  steady. 

8.  Mixing  is  completed  (the  flow  is  solved)  in  a  constant  area  (straight  tube)  duct. 

9.  The  surfaces  of  the  system  are  adiabatic. 

Two  parameters  that  will  be  varied  to  analyze  the  ejector  phenomena  are  entrained  mass 
flow  ratio  and  bypass  Mach  number.  Mass  flow  ratio  /i  is  defined  as 


ms 

M  =  — 
m* 

and  bypass  Mach  number  is  determined  from  the  pressure  ratio 


Pos 


(1) 


(2) 


Due  to  the  isentropic  relations  for  a  perfect  gas,  Pos  is  simply  related  to  Pi  by  specific 
heat  ratio  (constant)  and  the  Mach  number  of  the  secondary  flow.  Thus,  secondary'  flow 
Mach  number  (bypass  Mach  number),  is  the  only  variable  relating  Ps  to  Pos- 


1.2.3:  Free  Stre.^.m  Equations 

First  of  all,  one  must  know  the  stagnation  conditions  for  the  free  stream.  From  the 
definition  of  Mach  number  we  have 


VI  =  M\-sJy  R  T\ 


(3) 
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From  Table  2  we  know  M\  =  2.4  and  T\  =216.7  K;  since  the  fluid  is  air,  y  =  1 .4.  and  R 
287  J/(kg  ■  K),  so  this  gives 

V/  =  706.6  m/s 

According  to  the  stagnation  relations  for  temperature 

Tm  =  +  (4) 


Substituting  in  the  above  value  of  F/  from  Table  2  one  gets 
Toi  =  466.3  K 

The  isentropic  relations  also  give 

r  -1 

=  y 

Foi  \  Poi  / 
solving  for  Poi  one  gets 


Foi 


=  P\ 


r 

y-\ 


(5) 


(6) 


substituting  in  the  previously  calculated  values  gives 


=  169,500  Pa 
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1 .2.4:  Diffuser  Conditions 

For  the  preliminary  setup  of  the  spreadsheet,  the  conditions  of  station  2  and  3  are 
assumed  to  be  the  same.  Later,  a  straight  percentage  loss  in  total  pressure  will  be 
assumed  to  model  diffuser  losses.  But  for  the  initial  setup,  there  is  no  loss  in  total 
pressure  Po  between  station  1  and  5. 

1.2.5:  Ejector  (Second.ary  Flow)  Conditions 

For  the  ejector  (secondary  flow)  it  was  just  stated  that  its  stagnation  conditions  would  be 
assumed  to  be  the  same  as  the  freestream.  Therefore. 


Po5  =  Po,  (7) 

then  rewriting  (2) 

P5  =  PrP05  (8) 


solving  for  Mach  number  from  (4)  and  (5),  and  using  (2),  but  with  subscript  5  instead  of 
1 .  bypass  Mach  number  becomes 


.7-1 

2 

|r-  1 

(9) 


again  the  stagnation  relations  for  the  ejector  gives  {similar  to  (5),  except  with  subscript  5 
not  1 } 


75 


(10) 
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by  definition  of  Mach  number 

V5  =  R  Ts 

from  the  perfect  gas  relation 


Pi 


(11) 


(12) 


b\’  the  previously  stated  definition  of  p.  and  given  m4  is  already  fixed 

•  • 

mi  =  /u  mi  (13) 

finally  from  the  definition  of  mass  flow,  and  the  previous  equation 


Ai  =  // 


mi 

P  5  V5 


(14) 


1 .2.6:  Primary  Flow  Conditions 

For  the  primaiy  flow,  it  was  previously  stated  that  the  flow  would  be  completely 
expanded;  therefore,  it  was  expanded  to  match  its  static  pressure  with  that  of  the 
secondary'  flow;  so 

P4  =  PrP05  =  P5  (15) 

Given  that  Po4  is  fixed  from  the  Boeing  data,  and  using  the  just  calculated  P4.  M4  can  be 
calculated  from 
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.\/4  = 


r  r-i  1 

2 

y  .  1 

/- 1 

(16) 


gi\  en  that  T04  is  fixed  by  the  engine  data,  T4  can  be  calculated  from  the  isentropic  relation 
to  pressure 

/-I 

7-4  -  ^ 

primary  flow  velocity  is  easily  determined  from  the  definition  of  Mach  number 
V4  =  R  Ta  (18) 

density  comes  from  the  perfect  gas  relation 


PA  = 


Pa 
R  Ta 


(19) 


finally,  is  also  easily  determined  from  the  definition  of  mass  flow 


,-l4 


triA 

p  4  V4 


(20) 


1.2.7:  Flow  Mixing 

In  order  to  complete  the  analysis,  the  two  flows  must  mix.  It  was  assume  that  the  two 
flows  have  mixed  completely  and  perfectly  by  the  time  they  reach  station  six.  Beyond  the 
losses  incurred  for  mixing  itself,  no  additional  frictional  losses  were  accounted  for  in  the 
process.  Also,  the  distance  needed  to  complete  mixing  is  not  determined.  To  reiterate 
assumption  8.  one  item  to  note  is  that  the  mixing  is  accomplished  in  a  constant  area  duct. 
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which  has  an  area  given  as  the  sum  of  the  primaiv’  and  secondaiv'  flow  areas,  as 
previously  calculated,  so 

A6  =  A,  +  Aj  (21) 

1.2. 7.1;  Continuity 

To  solve  the  mixing  problem,  the  equations  of  continuity,  momentum  and  energy  are 
needed.^  Continuity'  for  a  perfect  gas  states  that 

•  •  •  • 

OT6  =  W4  +  ms  =  W4  (1  +  /j)  (22) 

1.2. 7.2:  Energy 

Next,  conservation  of  energy  states  that  the  energy  flow  into  the  mixing  chamber  equals 
the  energy  flow  out,  so 

Eb  =  E^  +  is  (23) 

In  general,  energy  flow  is  equal  to  the  following  according  to  the  definition  of  total 
enthalpy 

E  =  m  h(  (24) 

total  enthalpy  is  defined  as 

h,^h+'/by  =  CpT+'/2V-  (25) 


substituting  into  (24)  and  utilizing  the  definition  of  velocity  as  defined  by  Mach  number 
one  gets 
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•  • 

E  =  m 


CpT  +  -M-yRT 


(26) 


since  y  =  Cp/Cy 


E  =  m  CpT 


1  + 


M-  R 
2  C 


(27) 


plus  since  R  -  Cp  -  Cy,  and  remembering  the  definition  of  stagnation  temperature  similar 
to  equation  (4)  gives 


E  =  m  CpT 


1  +  ^ 
2 


=  m  CpTo 


(28) 


substituting  this  definition  of  energy  flow  (28)  into  equation  (23)  and  canceling  Cp  since 
it  is  assumed  to  be  constant  gives 

me  Toe  =  nn  To-i  +  ms  Tos  =  /W4  (7’o4  +  /uTos)  (29) 

this  can  easily  be  solved  for  Toe. 

1.2. 7.3:  Momentum 

This  leaves  the  momentum  equation  to  be  solved,  which  states  that  for  a  straight  tube 
solution 


•  •  • 

me  V6  +  Pe  Ae  =  {ot4  V4  +  Pa  Aa  +  ms  vs  +  Ps  As}  (30) 

Since  all  the  values  on  the  right  hand  side  of  the  equation  are  known,  for  simplicity  call 

•  •  • 

that  value  K.  (i.e.  K  =  /W4  V4  +  Pa  Aa  +  ms  vs  +  PsAs)  Mixed  mass  flow  me  is  given 

by  equation  (22),  and  Ae  is  given  by  (21),  but  Pe  and  ve  are  still  variables  in  question.  To 
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solve  this  dilemma,  first  rewrite  velocit>'  in  terms  of  Mach  number  and  using  the  perfect 
gas  law  along  with  the  definition  of  mass  flow  .  P(,  may  be  conveniently  eliminated. 

First  of  all,  since  mass  flow  \s  p  Ax  and  p  =  PI{R  T).  A  may  be  solved  for  and  is  equal  to 


m  RT 
Pv 


(31) 


substituting  in  the  definition  for  velocity  as  defined  by  Mach  number  (Equation  (3) 
without  subscript  1 ),  one  gets 


A  = 


m  RT 
PM  yjy  RT 


m  [¥Y 
PM\  y 


(32) 


Thus,  putting  (32)  back  into  (30),  but  with  subscript  6.  and  remembering  that  the  right 
hand  side  is  now  K,  one  gets 


•  Pe  me  R  Te 

me  V6  +  -  -  =  K 

PeMei  y 


(33) 


again  utilizing  the  definition  of  Mach  number  and  also  using  the  stagnation  temperature 
relation  for  Tt  (  similar  to  Equation  (4),  but  with  subscript  6) 


me 


Toe 


l  +  ^Me- 


Me 


1  R 


MeVy 


K 


(34) 


Equation  (34)  is  then  simply  a  function  of  the  mixing  plane  Mach  number  Me.  so  utilizing 
the  Newton  solver,  which  comes  with  the  EXCEL  spreadsheet,  a  solution  was  found. 
There  are  two  possible  solutions  to  the  mixing  plane  Mach  number.  These  are  the  so 
called  first  and  second  solutions.  The  first  solution  is  the  subsonic  solution,  and  the 
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second  solution  is  the  supersonic  one.  It  is  possible  to  get  only  one  solution  (i.e.,  it 
mechanically  chokes,  but  this  will  be  discussed  further  in  the  analysis  of  the  results).  For 
an  alternate  solution  to  this  problem,  where  M6  is  solved  for  explicitly,  see  Appendix  B."^ 
In  this  appendix,  an  analytic  solution  for  is  shown  and  the  conditions  under  which 
mechanical  choking  can  occur  are  summarized. 

1 .2.8:  Thrust  Generation  and  Flow  Expansion 

Now  that  the  flow  has  mixed,  it  must  be  “expanded”  back  to  ambient  pressure.  The  exact 
process  that  the  expansion/compression  cycle  would  take  is  dictated  by  the  end  of  mixing 
Mach  number.  There  are  three  possible  ways  to  bring  it  back  to  local  ambient  pressure. 
The  three  possibilities  are,  a  converging  duct,  a  diverging  duct,  and  finally,  a  converging 
then  diverging  duct.  One  of  these  three  will  allow  the  mixed  flow  to  “expand”  back  to 
the  current  ambient  pressure. 

Normally,  the  thrust  equation  for  an  open  cycle  engine  takes  the  following  form  (Which 
would  be  the  case  if  the  flow  were  not  “expanded”  through  a  system  of  ducts  and  was 
allowed  to  exit  at  the  mixing  cross  sectional  area). 

F  =  me  (\’6  -  Vr.)  +  Ae  {Pe  -  Pa>)  (35) 

However,  since  the  mixed  flow  will  be  expanded  back  to  ambient  pressure,  then  the 
pressure  at  the  "nozzle”  exit  will  be  equal  to  P^,  so  the  second  term  is  zero.  Recalling 
from  Figure  1.  that  the  nozzle  exit  plane  is  7,  thrust  then  becomes 


F  =  mi  (v7  -  \v.) 


(36) 
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Now  it  is  just  a  matter  of  determining  V7,  since  by  continuiw,  me  =  nv ,  and  =  v;  is 
already  known.  To  determine  v-.  we  need  to  know  Po-.  so  we  can  calculate  M?  and  then 
calculate  v;  from  there.  The  formulas  needed  to  do  this  are 


Po-  —  Poe  — 


P-  =  P. 


r 

=  P6  ^1  +  Me'Y  -  1 


A/’  = 


1 - 

1 

1  - 

2 

.1 

y-  1 

(37) 

(38) 

(39) 


and  since  the  process  is  an  isentropic  expansion,  stagnation  temperature  is  conserved,  so 
To7=Tn6  (40) 


V-  = 


Y  RToi  Ml' 


1  + 


r  -  1 


(41) 


Mi- 


Now,  having  calculated  v-,  the  thrust  produced  from  the  primary  flow/ejector 
combination  can  be  calculated  using  (36).  The  form  of  vy  in  Equation  (41)  is 
cumbersome  to  look  at,  but  it  is  implemented  that  way  in  the  spreadsheet  to  avoid  the 
calculation  of  Tj. 


1 .2.9;  Entropy  Production 

Knowing  there  are  two  possible  Mach  numbers  for  each  solution,  the  natural  question  is. 
"Which  one  will  be  the  one  that  will  naturally  occur?".  This  question  may  be  answered 
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by  looking  to  the  Second  Law  of  Thermodynamics. 

In  any  thermodynamic  process,  there  is  an  associated  entropy  change.  In  all  natural 
processes,  this  change  occurs  as  an  entropy  increase.  Thus,  if  the  process  would  have  a 
"negative"  change  in  entropy,  it  would  not  be  naturally  realizable.  So.  in  looking  at  the 
ejector  problem,  the  entropy  change  for  each  given  mixing  scenario  must  be  considered. 

From  the  well  known  thermodynamic  result,  for  an  adiabatic  system^ 

5  =  Cp  In  r  -  i?  In  /*  +  const  (42) 

one  can  conclude  that  the  total  entropy  for  an  open  system  is 

S  =  m  {Cp  In  T  -  R  \n  P  +  const)  (43) 

However,  since  we  are  assuming  a  thermally  and  calorically  perfect  gas.  the  constant  will 
be  dropped  assuming  the  same  thermodynamic  reference  for  all  gases. 

The  entropy  of  the  ejector  system  is  comprised  of  the  entropy  of  the  primary  flow  and  the 
secondary  flow  (stations  4  and  5  respectively).  Then  entropy  change  for  each  is  then 

Secondary  Flow 

AS2=  Si-Ss  Si  m5{Cp  In  Tb-R  In  Pb)-mi{Cp  In  Ti-R  In  Ps)  =  mi{Cp  ln—-R  In—) 

Ts  Pi 

Primary  Flow 

ASi  =  S'6-5'4  *  mi(Cp  In  Te-  R  In  Pb)-mi{Cp  In  Ti-  R  In  Pa)  =  mA{Cp  In— -  i?  In—) 

T  A  Pa 


So  the  total  entropy  change  is 


Referencing  it  to  the  primar\^  flow  by  dividing  through  by  niA  ■  R ,  one  gets 


\Sritt 

OT4  R 


A5i  + A52 
niA  R 


Cp  .  Tb  Cp  Tb  Pb  Pb 

=  — In —  +  //— In - In - /yin — 

R  Ta  R  Ts  Pa  Pi 


(45) 


Remembering  y  =  Q/C,-  and  R  =  Cp  -  Cv  ones  gets 


AiSIb/  Cp  .  Pb  Cp  I  Cv  .  Pb  .  Pb  ,  Pb 

-  =  - In — + // - In - In - am  — 

Cp-C.  Pa  RICv  Ps  Pa  Pi 


ASlnt 
nxA  R 


Y  Pb  Y  Pb  Pb  Pb 

- In —  +  a— — In - In— - //In  — 

Y-\  Pa  Y-'^  Ts  Pa  Pi 


But  since  the  primary  flow  is  fully  expanded,  and  P4  =  Pi 


r-1  ^4  y-i  r.  '  ’  Pa 


(46) 


(47) 


Therefore,  by  implementing  this  non-dimensionalized  total  entropy  check  into  the 
spreadsheet  analysis,  it  can  be  determined  what  solutions  will  be  physically  and 
thermodynamically  consistent  (i.e..  ASw/  >  0) 


CHAPTER  2:  INITIAL  ANALYSIS 


2.1:  SPREADSHEET  IMPLEMENTATION 

By  using  the  equations  from  the  previous  chapter,  a  spreadsheet  was  created  using 
Microsoft  EXCEL  5.0.  One  sheet  was  dedicated  to  each  value  of  Pr,  and  is  labeled  with 
the  appropriate  value  of  bypass  Mach  number.  On  each  sheet,  mass  flow  ratio  was  varied 
from  0.0  to  2.0  by  steps  of  0. 1 .  Then,  after  solving  the  momentum  equation,  and  solving 
the  expansion,  thrust  for  each  mass  flow  ratio  was  compared  to  the  baseline  {/u  =  0.0)  on 
each  sheet  and  was  expressed  as  a  percent  increase  or  decrease. 

Finally,  all  the  data  on  thrust  from  each  sheet  {Pr  =  0.9  -  0.4,  or  =  0.39  -  1.22).  was 
combined  on  one  sheet  plotting  bypass  Mach  number  versus  mass  flow  ratio  versus 
percent  thrust  increase.  To  determine  overall  thrust  augmentation,  all  values  were 
compared  to  value  of  thrust  calculated  for  mass  flow  ratio  equal  to  0.0.  This  is  simply  the 
thrust  of  the  primary  core  without  any  additional  ejector  flow. 

An  example  of  one  of  the  sheets  produced  by  EXCEL  is 


Table  3:  Worksheet  for  Pr  =  0.9003 
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One  thing  to  note  from  the  Worksheet  of  the  previous  page  is  the  column  labeled  Solver. 
This  is  the  column  containing  the  momentum  equation  to  be  solved.  The  equation  to  be 
solved  is  in  the  form  of  equation  (34)  from  Chapter  1.  The  Newton  solver  that  comes 
with  Microsoft  EXCEL  5.0  was  used  to  solve  for  a  solution.  One  thing  of  note  is  that 
seeded  properly  (initial  guess  value),  the  solver  can  converge  to  either  the  first  or  second 
solution.  In  order  to  get  a  solution,  and  initial  seed  guess  must  be  introduced  into  the 
solver,  and  it  is  from  this  seed  that  it  progresses  to  use  the  Newton  solver.  So.  in  the 
process  of  solving  the  problem,  careful  attention  was  paid  to  which  seed  value  was  used. 
Because,  for  example,  and  initial  guess  of  a  supersonic  end  of  mixing  Mach  number 
would  preclude  the  solver  from  being  able  to  converge  to  a  subsonic  root,  and  only 
Second  solution  answers  would  be  obtained. 

Another  item  of  interest  is  the  first  row  of  mass  flow  ratio.  The  0.0  mass  flow  ratio  is 
simply  the  situation  where  the  ejector  is  closed.  This  is  exactly  the  same  as  the  primaiy- 
core  acting  independently.  However,  as  a  consequence  of  this  spreadsheet  being  set  up  to 
handle  multiple  values  of  /u  and  Pr,  different  values  for  iV/g  will  show  up  as  the  solution 
to  the  fj.  —  0.0  ease,  for  various  values  of  Pr-  But  this  is  not  important,  since  the  mixing 
plane  for  the  zero  mass  flow  ratio  case  is  actually  non-existent.  The  “expansion”  that  it 
undergoes  is  a  “false”  one,  and  is  simply  a  consequence  of  the  way  the  table  is  set  up. 
The  important  thing  to  note  is  that  the  value  of  Mj  is  the  same  in  all  cases,  and  thus  the 
value  of  primary  core  thrust  remains  unchanged  while  varying  P^ 


2.2:  FIRST  VERSUS  SECOND  SOLUTION  COMPARISON 

When  looking  at  the  solution  to  a  supersonic  cruise  ejector  problem,  it  has  been  noted 
that  there  are  two  solutions  of  interest.  They  are  the  subsonic  mixing  solution  and  the 
supersonic  mixing  solution.  In  order  to  analyze  this  problem  completely,  both  solutions 
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must  be  compared.  They  will  be  compared  by  looking  at  entropy  change  and  percent 
thrust  increase. 

2.2.1:  First  Solution  D.ata 

Since  every'  entropy  change  calculated  would  be  excessively  cumbersome  to  present,  and 
not  entirely  informative,  selected  values  will  be  compared  from  each  solution.  The 
entropy  changes  that  will  be  looked  at  are  those  for  mass  flow  ratios  0,  1 .0  and  2.0  for 
values  of  P r  —  0.9,  0.7,  0.55  and  0.4.  The  values  of  correspond  to  bypass  Mach 
numbers  of  0.j9,  0.73,  0.97  and  1.22  respectively.  The  above  choices  of  points  to 
evaluate  are,  of  course,  completely  arbitrary',  but  give  12  points  of  comparison  from  a 
possible  240  (the  number  of  data  points  in  each  thrust  table  used  to  create  the 
forthcoming  3D  graphs). 

With  regard  to  the  percent  thrust  increase,  all  240  points  will  be  shown  in  both  the 
subsonic  and  the  supersonic  case. 

The  values  for  the  entropy  change,  for  the  first  solution  to  the  mixing  problem  at  the 
previously  indicated  points  are: 

Table  4:  First  Solution  Non-dimensional  Entropy 
Changes 


Values  of 

ASinl 

Pr  =  0.9003 

Pr  =  0.70 

Pr  =  0.55 

II 

O 

• 

/W4  R 

p=o:o 

3.28E'" 

-2.8E‘‘^ 

5.55E-'^ 

-5.6E‘'^ 

M=  1.0 

0.924851 

0.856494 

0.862277 

0.947219 

11 

b 

1.359347 

1.281898 

1.216471 

1.361691 
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The  values  of  percent  thrust  increase  were  all  calculated  by  referencing  the  =  0.0,  = 

0.9003  case.  For  the  first  solution,  the  values  of  percent  increase  in  thrust  are  shown  in 
the  following  table  and  graphically  displayed  in  the  subsequent  figure. 


Table  5;  First  Solution  Thrust  Augmentation 


Thrust  when  expanded  to  ambient  Pressure  - 
%  Thrust  increase  over  baseline,  -107  kN 


vJIPr^ 

0.90 

0.85 

0.80 

0.75 

0,70 

0.65 

0.60 

0.55 

0.5328 

0.50 

0.45 

0.40 

0.39 

0.49 

0.57 

0.65 

0.73 

0.81 

0.89 

0.97 

1.0000 

1.05 

1.13 

1.22 

0.0 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.1 

0.88 

0.97 

1.00 

0.97 

0-88 

0.72 

0.48 

0.14 

0.00 

-0.31 

-0.90 

-1.67 

0.2 

1.71 

1.92 

2.01 

2.04 

2.00 

1.90 

1.70 

1.40 

1.27 

0.98 

0.40 

-0.37 

0.3 

2.47 

2.78 

2.92 

2.99 

2.99 

2.93 

2.77 

2.51 

2.39 

2.11 

1.54 

0.77 

0.4 

3.19 

3.56 

3.75 

3.84 

3.88 

3.84 

3.72 

3.49 

3.37 

3.11 

2.56 

1.78 

0.5 

3.85 

4.28 

4.50 

4.62 

4.66 

4.65 

4.56 

4.36 

4.26 

4.00 

3.46 

2.69 

0.6 

4.46 

4.94 

5.19 

5.32 

5.37 

5.38 

5.32 

5.14 

5.04 

4.81 

4.28 

3.50 

0.7 

5.03 

5.55 

5.81 

5.96 

6.05 

6.18 

5.99 

5.84 

5.76 

5.53 

5.01 

4.24 

0.8 

5.56 

6.11 

6.39 

6.54 

6.71 

6.91 

6.77 

6.48 

6.40 

6.19 

5.68 

4.90 

0.9 

6.05 

6.63 

6.93 

7.09 

7.31 

7.58 

7.49 

7.06 

6.98 

6.78 

6.29 

5.51 

1.0 

6.51 

7.11 

7.42 

7.59 

7.85 

8.20 

8.15 

7.74 

7.52 

7.33 

6.85 

6.07 

1.1 

6.94 

7.56 

7.88 

8.05 

8.35 

8.76 

8.76 

8.38 

8.16 

7.83 

7.36 

6.58 

1.2 

7.34 

7.98 

8.31 

8.48 

8.80 

9.29 

9.33 

8.97 

8.75 

8.28 

7.83 

7.04 

1.3 

7.72 

8.37 

8.71 

8.89 

9.21 

9.77 

9.85 

9.52 

9.30 

8.77 

8.26 

7.48 

1.4 

8.07 

8.74 

9.09 

9.27 

9.59 

10.21 

10.34 

10.03 

9.82 

9.28 

8.66 

7.88 

1.5 

8.41 

9.08 

9.44 

9.62 

9.93 

10.63 

10.80 

10.51 

10.30 

9.76 

9.04 

.  8.25 

1.6 

8.72 

9.41 

9.77 

9.95 

10.25 

11.01 

11.23 

10.96 

10.75 

10.21 

9.39 

8.60 

1.7 

9.02 

9.71 

10.08 

10.26 

10.54 

11.37 

11.63 

11.38 

11.18 

10.64 

9.71 

8.92 

1.8 

9.30 

10.00 

10.37 

10.56 

10.81 

11.70 

12.00 

11.77 

11.58 

11.04 

10.01 

9.22 

1.9 

9.56 

10.27 

10.65 

10.84 

11.05 

12.02 

12.36 

12.15 

11.95 

11.41 

10.29 

9.50 

2.0 

9.81 

10.53 

10.91 

11.10 

11.28 

12.31 

12.69 

12.50 

12.31 

11.76 

10.56 

9.77 

Two  things  to  note  about  the  previous  table  are,  first,  the  boxed  values  indicate  that  the 
momentum  equation  could  either  not  be  solved,  or  choked  flow  occurred.  To  determine 
if  choking  resulted  or  if  no  solution  occurred  was  a  matter  of  looking  at  the  column 
labeled  "Solver’  as  was  shown  in  Table  3.  If  the  end  of  mixing  Mach  number  was 
calculated  to  be  1  or  close  to  it,  and  the  value  of  the  solver  was  zero,  then  only  one 
solution  was  found  (the  choked  solution).  But,  for  all  of  the  cases  ran,  that  combination 
of  factors  rarely  occurred  for  a  Mach  number  exactly  equal  to  1  (However  the  end  of 
mixing  Mach  number  converged  to  1  even  if  the  solver  was  non-zero). 
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So.  for  all  of  the  boxed  values  above,  these  values  of  thrust  augmentation  are  not 
possible,  since  the  momentum  equation  could  not  be  solved  for  a  given  Mach  number  in 
the  mixing  plane.  Thus,  any  value  in  a  table  of  thrust  augmentation  boxed  with  a  single 
line,  is  not  considered  a  viable  point  where  thrust  augmentation  could  occur. 

The  second  thing  to  note  about  the  previous  table  is,  the  values  of  mass  flow  ratio  are  on 
the  vertical  axis  of  the  graph  and  the  horizontal  axis  has  pressure  recoveiv'  Pr,  with  bypass 
Mach  number  underneath. 


%lncreas8  Thrust,  expanded  to  ambient 


Bypass 

Mach 

Number 


□ 

12.00-15.00 

■ 

9.00-12.00 

Mass 

□ 

6.00-9.00 

How 

Ratio 

□ 

3.00-6.00 

■ 

0.00-3.00 

■ 

-3.00-0.00 

Figure  2;  First  Solution  Graph  of  %  Thrust 
Augmentation 


2.2.2:  Second  Solution  Dat.a 

For  the  second  solution,  there  are  some  noticeable  changes  from  the  first  solution.  First 
of  all.  it  is  possible  to  violate  the  Second  Law  of  Thermodynamics,  and  secondly,  the 
thrust  augmentation  achieved  is  much  larger,  for  similar  conditions. 

The  entropy  values  corresponding  to  the  same  points  as  was  presented  in  the  first  solution 


are 
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Table  6:  Second  Solution  Non-dimensional 
Entropy  Changes 


Values  of 

A*Sfo/ 

Pr  =  0.9003 

?,  =  0.70 

Pr  =  0.55 

Pr  =  0.4 

nu  R 
//=0.0 
//=  1.0 
/^=2.0 

4.44E-'^ 

0.301295 

-0.64536 

-2.8E-'* 

0.856494 

1.281898 

5.55E'‘^ 

0.862277 

1.216471 

-5.6E''^ 

0.904364 

1.341958 

The  table  of  thrust  augmentation  for  the  second  solution  shows  a  definite  trend  towards 
larger  thrust  increases.  Yet  also  note  that  there  are  double  outlined  boxes  which  indicate 
the  solution  would  violate  the  Second  Law,  and  as  before  the  single  outline  indicates 
solutions  which  do  not  solve  the  momentum  equation.  The  corresponding  table  and 
figure  for  the  second  solution  follow: 


Table  7:  Second  Solution  Thrust  Augmentation 


Thrust  when  expanded  to  ambient  Pressure 
%  Thrust  increase  over  baseline,  -107  kN 


Ml! 

0.90 

0.85 

0.80 

0.75 

0.70 

0.65 

0.60 

0.55 

0.5328 

0.50 

0.45 

0.40 

0,39 

0.49 

0.57 

0.65 

0.73 

0.81 

0.89 

0.97 

1.0000 

1.05 

1.13 

1.22 

0.0 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.1 

1.26 

1.15 

1.14 

1.15 

1.15 

1.14 

1.13 

1.11 

1.10 

1.09 

1.05 

1.01 

0.2 

2.69 

2.27 

2.18 

2.17 

2.17 

2.16 

2.14 

2.11 

2.10 

2.06 

2.00 

1.93 

0.3 

4.28 

3.35 

3.13 

3.09 

3.09 

3.08 

3.06 

3.01 

2.99 

2.95 

2.87 

2.76 

0.4 

5.98 

4.41 

4.01 

3.93 

3.93 

3.92 

3.89 

3.84 

3.81 

3.76 

3.65 

3.51 

0.5 

7.79 

5.45 

4.83 

4.69 

4.68 

4.68 

4.65 

4.59 

4.56 

4.49 

4.37 

4.21 

0.6 

9.70 

6.49 

5.61 

5.39 

5.38 

5.38 

5.35 

5.27 

5.24 

5.17 

5.02 

4.84 

0.7 

11.69 

7.51 

6.34 

6.04 

6.05 

6.18 

5.99 

5.91 

5.87 

5.79 

5.63 

5.42 

0.8 

13.77 

8.54 

7.04 

.  6.64 

6.71 

6.91 

6.77 

6.50 

6.45 

6.36 

6.19 

5.96 

0.9 

15.92 

9.56 

7.71 

7.20 

7.31 

7.58 

7.49 

7.06 

7.00 

6.89 

6.70 

6.45 

1.0 

18.14 

10.58 

8.36 

7.72 

7.85 

8.20 

8.15 

7.74 

7.52 

7.39 

7.18 

6.91 

1.1 

20.43 

11.60 

8.98 

8.21 

8.35 

8.76 

8.76 

8.38 

8.16 

7.85 

7.63 

7.34 

1.2 

22.78 

12.63 

9.59 

8.67 

8.80 

9.29 

9.33 

8.97 

8.75 

8.29 

8.04 

7.74 

1.3 

25.18 

13.66 

10.18 

9.11 

9.21 

9.77 

9.85 

9.52 

9.30 

8.77 

8.43 

8.11 

1.4 

27.64 

14.69 

10.76 

9.53 

9.59 

10.21 

10.34 

10.03 

9.82 

9.28 

8.79 

8.46 

1.5 

30.14 

15.73 

11.33 

9.93 

9.93 

10.63 

10.80 

10.51 

10.30 

9.76 

9.13 

8.79 

1.6 

32.69 

16.77 

11.89 

10.31 

10.25 

11.01 

11.23 

10.96 

10.75 

10.21 

9.45 

9.09 

1.7 

35.29 

17.82 

12.44 

10.68 

10.54 

11.37 

11.63 

11.38 

11.18 

10.64 

9.76 

9.38 

1.8 

37.93 

18.88 

12.98 

11.04 

10.81 

11.70 

12.00 

11.77 

11.58 

11.04 

10.04 

9.65 

1.9 

40.61 

19.94 

13.52 

11.38 

11.05 

12.02 

12.36 

12.15 

11.95 

11.41 

10.31 

9.90 

2.0 

43.32 

21.00 

14.05 

11.71 

11.28i 

12.31 

12.69 

12.50 

12.31 

11.76! 

10.57 

10.14 

%lncrease  Thrust,  expanded  to  ambient 


Bypass 

Mach 

Number 


■  40.00-45.00 

□  35.00-40.00 

■  30.00-35.00 

□  25.00-30.00 

■  20.00-25.00 

□  15.00-20.00 

□  10.00-15.00 

■  5.00-10.00 

■  0.00-5.00 


Figure  3:  Second  Solution  Graph  of  %  Thrust 
Augmentation 


Looking  at  the  graph  above,  one  would  be  tempted  to  say  that  low  bypass  Mach  numbers 
would  be  the  best  solution,  due  to  the  large  spike  in  thrust  augmentation.  However,  this 
would  be  a  false  conclusion,  due  to  the  violation  of  the  Second  Law  that  is  indicated  by 
the  double  outlined  values  in  Table  7.  Yet.  the  second  solution  still  seems  promising, 
since  it  has  higher  thrust  augmentation  then  the  first  solution. 


But.  without  an  exact  comparison  of  mixing  lengths  required,  nothing  can  be  said  about 
which  solution  is  better.  Assuming  a  similar  rate  of  mixing,  the  supersonic  solution 
would  take  a  longer  distance  to  mix  due  to  its  greater  flow  velocity.  Since  nothing  can  be 
said  about  which  solution  provides  more  realizable  thrust  augmentation,  without  a  valid 
mixing  model,  both  results  will  be  presented. 


2.3:  CONCLUSIONS 

So.  by  looking  at  the  first  and  second  solution,  one  would  want  to  say  that  the  second 


solution  is  unequivocally  better,  by  virtue  of  thrust  augmentation  alone.  However,  no 
analysis  was  made  of  the  actual  distance  required  to  mix,  so  the  second  solution's 
apparent  benefits,  could  not  be  compared  to  the  possible  drawbacks  of  a  longer  mixing 
length.  The  thrust  augmentation  increase  may  be  offset  by  the  associated  weight  penalty 
of  a  longer  mixing  length.  There  is  room  for  further  analysis  in  this  area,  and  if  this  is  to 
be  done  a  model  for  mixing  must  be  made. 


CHAPTERS:  DIFFUSER  LOSSES 


After  looking  at  the  simple  case  of  the  ejector  in  cruise  configuration  for  a  subsonic  and 
supersonic  mixing  solution,  it  is  natural  to  extend  this  by  attempting  to  account  for 
diffuser  losses.  It  was  assumed  that  the  entrained  ejector  air  was  brought  in  through  the 
diffuser  without  any  loss  in  total  pressure.  For  the  case  without  losses,  this  was  all  right, 
but  assuming  the  ejector  diffuser  is  imperfect  and  may  have  some  losses  associated  with 
it.  losses  must  be  introduced.  This  was  done  by  simply  taking  a  percentage  loss  in 
stagnation  pressure  from  the  ambient  flow. 


3.1:  EFFECT  OF  DIFFUSER  LOSS 

This  subsequent  decrease  in  ejector  total  pressure,  affects  its  static  pressure  too.  which 
therefore  affects  the  full  expansion  of  the  primary  flow.  Losses  were  taken  in  5% 
increments  from  5%  to  25%.  The  affects  of  these  losses  are  shown  in  the  following 
tables  and  charts. 

3.1.1:  5%  Diffuser  Loss  (1st  AND  2nd  Solution) 

This  section  and  the  following  sections  will  have  information  on  both  the  first  and  the 
second  solution.  Tables  of  entropy  changes,  percent  thrust  increase  may  be  found  in 
Appendix  C;  however,  the  graphs  of  the  thrust  augmentation  will  be  included  in  each 


section. 
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%increase  Thrust,  expanded  to  ambient 


Bypass 

Mach 


Number 


■  9.00-12.00 

□  6.00-9.00 

□  3.00-6.00 

■  0.00-3.00 
■ -3.00-0.00 


Figure  4:  Graph  of  %  Thrust,  5%  Loss  (1st 
Solution) 


%lncrease  Thrust,  expanded  to  ambient 


Bypass 


Mach 

Number 


■  40.00-45.00 

□  35.00-40.00 

■  30.00-35.00 

Mass 

□  25.00-30.00 

Row 

■  20.00-25.00 

Ratio 

□  15.00-20.00 

□  10.00-15.00 

■  5.00-10.00 

■  0.00-5.00 

Figure  5:  Graph  of  %  Thrust,  5%  Loss  (2nd 
Solution) 


ro 
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.1.2:  10%  Diffuser  Loss  (1st  and  2nd  Solution) 


%lncrease  Thrust,  expanded  to  ambient 


Bypass 

Mach 

Number 


Mass 

Row 

Ratio 


□  6.00-9.00 

□  3.00-6.00 
■  0.00-3.00 
a -3.00-0.00 


Figure  6:  Graph  of  %  Thrust.  10%  Loss  (1st 
Solution) 

%lncrease  Thrust,  expanded  to  ambient 


Bypass 

Mach 

Number 


■  36.00-42.00 

□  30.00-36.00 

■  24.00-30.00 
018.00-24.00 

□  12.00-18.00 

■  6.00-12.00 
■  0.00-6.00 


Figure  7:  Graph  of  %  Thrust.  10%  Loss  (2nd 
Solution) 
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3.1.3;  15%  Diffuser  Loss  (1st,\nd  2nd  Solution) 


“/cincrease  Thrust,  expanded  to  ambient 


Bypass 


Mach 


Mass 

Row 

Ratio 


Number 


■  4.00-6.00 
n  2.00-4.00 
□  0.00-2.00 
■ -2.00-0.00 
a  -4.00-2.00 


Figure  8:  Graph  of  %  Thrust,  15%  Loss  (1st 
Solution) 


%lncrease  Thrust,  expanded  to  ambient 


Bypass 

Mach 

Number 


Mass 

Row 

Ratio 


□  30.00-40.00 

□  20.00-30.00 
■  10.00-20.00 
■  0.00-10.00 


Figure  9:  Graph  of  %  Thrust,  15%  Loss  (2nd 
Solution) 
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3.1.4:  20%  Diffuser  Loss  (1st  .\nd  2nd  Solution) 


Voincrease  Thrust,  expanded  to  ambient 


Bypass 

Mach 


■  2.00-3.00 
01.00-2.00 

■  0.00-1.00 
□ -1.00-0.00 
□ -2.00-1. 00 
■ -3.00-2.00 
0-4.00-3.00 


Number 


Figure  10:  Graph  of  %  Thrust,  20%  Loss  (1st 
Solution) 


Vclncrease  Thrust,  expanded  to  ambient 


Bypass 

Mach 

Number 


□  35.00-40.00 

■  30.00-35.00 
0  25.00-30.00 

■  20.00-25.00 

□  15.00-20.00 

□  10.00-15.00 

■  5.00-10.00 

■  0.00-5.00 


Figure  1 1 :  Graph  of  %  Thrust,  20%  Loss  (2nd 
Solution) 
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3.1.5:  25%  Diffuser  Loss  (1st  AND  2nd  Solution) 


%lncrease  Thrust,  expanded  to  ambient 


Bypass 

Mach 

Number 


■  1.00-2.00 

□  0.00-1.00 

Mass 

■ -1.00-0.00 

Row 

□ -2.00-1.00 

Ratio 

□ -3.00-2.00 

■ -4.00-3.00 

H -5  00-4.00 

Figure  12:  Graph  of  %  Thrust,  25%  Loss  (1st 
Solution) 


%lncrease  Thrust,  expanded  to  ambient 


Bypass 

Mach 

Number 


Mass 

Row 

Ratio 


■  35.00-40.00 

□  30.00-35.00 

■  25.00-30.00 

□  20.00-25.00 

■  15.00-20.00 

□  10.00-15.00 

□  5.00-10.00 

■  0.00-5.00 
■ -5.00-0.00 


Figure  13:  Graph  of  %  Thrust,  25%  Loss  (2nd 
Solution) 
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3.2:  COMMENTS  ON  DIFFUSER  LOSS 

Looking  at  the  previous  subsections  one  can  see  an  overriding  trend  of  thrust 
augmentation  decrease.  As  the  diffuser  becomes  less  efficient  at  recovering  total  pressure 
from  the  ambient  flow,  the  thrust  of  the  ejector  combination  begins  to  decline. 

Once  again  the  second  solution  has  a  larger  thrust  augmentation,  but  once  again  certain 
parts  are  defeated  by  the  Second  Law  of  Thermodynamics.  While  nothing  can  be  said 
about  which  solution  is  easier  to  achieve  practically  (the  first  or  the  second),  they  both 
show  a  positive  non-dimensional  entropy  change  over  a  large  range  of  values.  So,  both 
solutions  are  of  interest,  and  should  be  examined. 

Also,  the  tables  of  thrust  augmentation  still  show  a  problem  with  the  momentum  equation 
not  being  satisfied.  However,  one  interesting  trend  is  the  disappearance  of  this 
phenomena  with  increased  diffuser  loss.  The  flow  becomes  closer  to  being  choked  (or 
sonic),  and  this  can  be  explained  physically,  by  the  following.  The  choking  occurs,  since 
the  primary  stream  is  trying  to  expand  to  the  secondary’  stream  static  pressure.  In  doing 
so,  there  are  certain  configurations  of  mass  flow  ratio  and  bypass  Mach  number  (the 
single  boxed  values)  where  the  primary  flow  continues  to  expand  in  the  ejector,  but  leads 
to  an  effective  constriction  of  the  area  the  secondary  flow  must  pass  through.  This 
restricted  area  acts  as  a  virtual  throat  for  the  secondary'  flow  to  pass  through.  Regardless 
of  whether  the  bypass  (secondary)  flow  is  initially  subsonic  or  supersonic,  the  restriction 
in  area  acts  to  bring  the  flow  back  to  sonic.  A  similar  process  is  simultaneously 
happening  with  the  primary’  flow,  so  it  goes  sonic  too.  With  increasing  losses  in 
secondary  stream  stagnation  pressure,  there  is  also  an  associated  loss  in  static  pressure, 
since  the  bypass  Mach  number  remains  the  same.  This  more  readily  allows  for  a  solution 
of  the  momentum  equation  that  is  not  choked. 


CHAPTER  4:  BLEED  LOSSES 


4. 1 :  TREATMENT  OF  BLEED  LOSS 

After  doing  the  analysis  given  so  far,  the  results  were  presented  to  the  Boeing  Company. 
During  the  meeting  with  Boeing,  Mr.  Jan  Syberg  noted  that  there  was  no  account  for 
bleed  losses.  This  led  the  author  to  consider  bleed  losses  in  his  analysis,  and  to  model 
something  similar  to  what  Boeing  said  was  feasible,  a  six  percent  diffuser  loss  along  with 
a  five  percent  bleed  from  the  ejector  was  what  was  used  for  the  model. 


An  attempt  to  model  bleed  loss  was  done  by  treating  it  as  a  throttling  process.  It  was 
assumed  that  the  analysis  done  previously  would  be  accurate  for  the  flow  that  was  not 
bled  from  ejector’s  diffuser.  So,  since  the  previous  analysis  valid,  the  bleed  flow  could 
be  handled  separately.  The  bleed  flow  was  calculated  as 


OTb 


. 

ms 


ms 

J 


ms  (0.05263) 


(48) 


This  quantity  was  fixed  by  continuity,  but  for  the  process  of  throttling,  the  bleed  flow  was 
considered  to  be  adiabatic,  isentropically  expanded  and  then  shocked  back  to  a  lower 
pressure.  For  one  case  it  was  shocked  back  to  atmospheric  pressure,  and  for  the  second 
case,  it  was  throttled  to  lose  one-half  of  its  stagnation  pressure. 

The  figure  of  the  next  page  shows  how  the  process  for  the  first  case  was  modeled: 


/C 

Ab,^  ^ 


Isentropic  Expansion 


-  Normal  Shock 
Ah2 _ 


Ejector 


Engine 


Figure  14:  Throttling  Process  Diagram  for 
Throttle  to  Ambient 

The  value  for  was  calculated  by  the  following  equation. 


^  mb  mb 

Ab\  =  -  =  - 

/7b  Vb  p\V\ 


(49) 


This  assumes  that  the  inlet  velocity  for  the  bleed  flow  is  exactly  the  same  as  the 
translational  velocity  of  the  aircraft. 


For  the  second  case,  an  attempt  was  made  to  bleed  more  intelligently.  This  led  to  the 
result  of  losing  1/2  Po-  To  do  this  the  bleed  flow  was  bled  just  prior  to  the  point  in  the 
diffuser  where  the  flow  was  sonic,  i.e.  .V/  =  1 .  From  there  it  would  be  shocked  and  then 
expanded  to  ambient  pressure.  Across  this  shock,  from  the  Mach  2.4  flow,  the  ratio  of 
stagnation  pressure  after  the  shock  (Poj)-  to  that  before  the  shock  (Poi)  is  P n2^P oi  ~  0.5.^ 
So  the  choice  was  made  to  use  1/2  Po  as  the  model  for  bleed  in  the  second  case.  A 
diagram  of  the  process  of  the  second  case  is  the  following: 
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-  Analogy  to 

inlet  diffuser 
shock 

Expand  isentropically  to  ambient  pressure 


■  Bleed  from  throat  of  ejector  (outside  diagram  is  a  model) 


Engine 


Figure  15:  Throttling  Process  Diagram  for 
Throttle  to  1/2  Po 

For  purposes  of  the  following  derivations  both  equations  (48)  and  (49),  that  were  given 
previously,  hold. 

4.1.1:  Throttle  to  Atmospheric  Pressure 

For  the  first  case,  it  was  desired  that  the  pressure  after  the  shock  would  equal  atmospheric 
pressure.  So,  to  do  this  the  following  relations  and  variable  definitions  were  used.  Po  is 
the  stagnation  pressure  of  the  flow,  Pe  is  the  static  pressure  at  the  exit  of  the  throttle  prior 
to  shock,  and  P2  is  the  static  pressure  of  the  flow  after  the  shock. 

P.  =  + 
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Ml 


Me-  + 


2 


IMc- 


r 


r-\ 


-1 


Pz  \+rMe- 
Pe  ~  \  +  Y  Mz- 

Combining  these  gives: 

Pz  _  \  +  y  Me- 

2  ^ 

Me-  + - 

r-1 

2M'  -^-1 

V  y-\  J 


Po  \  l  +  - 

2  J 


y-i 


1  +  7 


(50) 


This  holds  for  any  throttling  process.  However,  for  the  first  case,  the  static  pressure  after 
the  shock  was  desired  to  be  atmospheric  pressure.  This  is  to  say  that  P2  =  Pi.  Another 
relation  that  will  simplify  the  equation  is  that  in  air  (where  y  =  1 .4),  Po  =  Vo  0684*f’/-  (Cite 
the  Emmons  Gas  dynamics  tables)  Substituting  these  back  into  (50)  and  performing 
some  algebra  gives: 


Pz 


1 

00684 


Pz 


Y 

7-1 


2M--  -^-1  +  7 
7-1 


2  ^ 

y-\j 


Canceling  the  P2,  noting  that  for  air,  =  .2.  7.,.i=3.5.  and  %.[=  5:  and  finally  replacing 
Me  with  M. 
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0.0684  (l  +  1.4A/^)(7A/--l) 

(l  +  0.2M^P  ”7M- -1  +  1.4  (M-+5) 

Solving  numerically  for  M  gives  =  4.714162  «  4.7.  Now  that  the  Mach  number  to 
which  the  flow  must  be  isentropically  expanded  is  known,  the  thrust  decrease  may  be 
calculated  from  the  following  formulas. 

Floss  -  Fml  +  FpL 


Fml  =  rnh  (v2  -  vi) 

F PL  =  {P 2  -Pi)  A b2 

where  F\/i  is  the  thrust  loss  due  to  momentum  decrease  and  Fpi  is  the  thrust  '‘loss"  due  to 
pressure  differences. 

For  this  first  case,  since  P:  =  Pi,  the  thrust  “loss”  due  to  pressure  is  zero,  and  this  only 
leaves  the  momentum  loss  term.  To  calculate  this  term,  the  value  of  must  be 
calculated.  To  calculate  this  velocity,  one  must  note  that  the  speed  of  sound  used  to 
calculate  is  different  from  that  of  the  one  used  to  calculate  v/.  This  is  due  to  the 

temperature  change  that  occurred  in  the  flow  while  it  was  being  isentropically  and 

) 

adiabatically  expanded  to  M=  4.7. 

Even  though  Ab2  will  not  be  used  in  this  first  case,  the  process  of  its  calculation  will  be 
covered,  since  it  will  actually  be  used  in  the  second  case.  The  normal  shock  relations  for 
M  =  4.7  that  will  be  used  to  calculate  Ab2  and  are  as  follows 

19.58  @M=  4.7 


A2.4lA*  =  2.m\  @M=2.4 
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TITo  =  0.46468  @  M  =  2.4 

r/ro  =  0.18457  @M=  4.7 

T:ITe  =  5.2334  (across  a  M=  4.7  normal  shock) 

Using  these  to  calculate  Ab2  gives 

Am  a*  ,  19.58  1  ,  , 

Jf,'  = - Ahi  = - Am  =  8.2602  .4/>i 

A*  A2  4  1  2.4031 

Now  to  calculate  v?,  the  new  static  temperature  must  be  calculated.  This  can  then  be 
substituted  into  the  following  equation  to  give  vj.  This  equation  can  also  be  used  to 
calculate  v/. 

V  =  Myir  R  T 

For  the  free  stream  velocity  one  gets  the  following 

VI  =  2.4Vl.4-(287  J/kg-K)-(216.7  K)  =  708.183  m/s 

Using  the  shock  relations  detailed  above,  the  static  temperature  after  the  shock  is 

T2  =  5.2334  T,  =  5.2334  '  T, '  ToilT, '  TJToe  =  5.2334-216.7-'/o.46468-°  =  450.4535  K 

So  in  this  new  flow  0.41992,  via  the  normal  shock  relation  for  M=  4.7,  and  using  the 
definition  of  velocity  above  one  gets 

V2  =  0.41 992^1.4 -(287) -(450.4535)  =  178.647  m/s 

This  then  gives  momentum  thrust  loss  as 

Fml  =  mh  {v2-Vi)  =  mh  (178.647  -  708.183)  =  m/.  (-529.5359  m/ s) 
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Now  with  these  new  values  calculated,  five  new  columns  were  added  to  the  previously 

used  spreadsheets.  These  columns  were  mb ,  Abi,  Ab2,  Fml  and  Fpl.  Adding  the  values  of 
momentum  thrust  and  pressure  thrust  into  the  previously  calculated  trust  gives  new  values 
for  overall  thrust  and  leads  to  a  change  in  the  percent  thrust  increase. 


4. 1 .2:  Throttle  to  lose  half  St.agnation  Pressure 

Due  to  the  reasons  stated  before,  the  throttle  to  1/2  Po  was  chosen  as  the  second  case. 
But.  in  order  to  evaluate  this  case,  a  slight  variation  on  the  equations  just  derived  must  be 
calculated.  While  equation  (50)  above  would  hold  for  this  case,  it  is  not  in  a  convenient 
form  to  relate  stagnation  pressures.  Therefore  the  following  derivation  was  used. 

P20  P20  P\  ,  P20  P20  P2 

- = - where - = - 

P\0  P\  P\0  Pi  P2  Pi 


El 

Pi 


\  +  y  Mi- 
1  +  7  Af2- 


r-i 


Since  we  want  P20  to  be  ’/2  Pio,  then 


1  _  P20  _  P2a  Pi  P2 

2  PiO  P2  PiO  Pi 


r 

r-\ 


\  +  r  Ml 

1  +  7  M2- 


Mi  H - 

7-1 


2M\- 

7-1 


since  M2 


,  substitution  into  the  above  with  y  =  1 .4  gives 
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2 


1  +  0.2 


M\-  +  5^ 
1M\-  -\) 


(1  +  0.2  Mr)"'' 


1  +  1.4  M\- 


1  +  1.4 


M\-  +5 
7Mr  -1 


After  performing  some  algebra,  and  letting  Mi  =  M,  one  gets 


2 


fl  + 


1  +  0.2.V/-^ 

7M-  -1 


1  +  0.2A/- 


V 


y 


(7.U'  -l)(l  +  1.4  M-) 
1 M-  -  \  +  \A[m-  +5) 


2 


|(l+0.2M')‘'  +(7.V/'  -l)'’j'' 


9.8  M'*  +5.6  A/-  -1 
8.4  M-  +6 


Numerically  solving  for  M.  one  finds  that  an  exit  Mach  number  of  .V/=  2.4975  =  2.5  prior 
to  the  shock  gives  a  decrease  in  stagnation  pressure  of  'A.  Thus,  the  throttling  process 
involves  an  expansion  from  M=  2.4  to  Af  =  2.5. 


Expanding  the  flow  to  Mach  2.5  creates  a  new  area  ratio  as  well  as  a  new  static 
temperature  and  therefore  speed  of  sound.  The  important  information  from  the  normal 
shock  tables  that  will  be  needed  to  calculate  Abi  and  T:  are 


A2.5IA*  =  2.61,61  @M=1.5 

T/To  =  0.44444  @M=  2.5 

TA/Te  =  2.1375  (across  a  M=  2.5  normal  shock) 


using  these  and  the  information  from  the  previous  section. 


Ah2  = 


A2.5  A* 
A  *  A2.4 


A/jI  = 


2.6367 


1 


Abi  =  1.0972  Am 


1  2.4031 
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T:  =  2.1375  T,  =  2.1375-  Tj  ’  ToilTi '  TJToe  =  2.1375-216.7-'/o.46468'°‘^'^'‘/i  =  443.0209  K 

However,  as  before  there  is  a  decrease  in  Mach  number  across  a  normal  shock  and  for  M 
=  2.5,  the  Mach  number  decreases  to  M=  0.5 1299  and  plugging  these  in  gives 

V2  =  0.5 1299V1.4- (287) -(443.0209)  =207.261  m/s 

This  then  gives  momentum  thrust  loss  as 

Fml  =  mh  (v:  -  vi)  =  mh  (207.261-708.183)  =  (-491.749  m/s) 


But.  one  more  piece  of  information  is  needed  to  calculate  Fpl-  The  pressure  on  the  other 
side  of  the  shock  is  needed.  This  can  be  calculated  in  terms  of  the  initial  inlet  static 
pressure  since  stagnation  pressure  is  conserved  up  to  the  point  of  the  shock,  and  then 
there  is  a  relation  betw'een  the  static  pressures  before  and  after  the  shock.  For  the 
isentropic  expansion  to  Me  =  2.5 


y 

Pe  =  Po{\  +  ^^  Mc-y''  ^Po  (0.05853)  = 


Pi  (0.05853) 
0.0684 


=  0.8557  P. 


but  across  a  Mach  2.5  normal  shock 
PiJPe  =  7.125  so  P.  =  7.125  Pe  =  6.0969  P/ 

now,  substituting  all  of  these  into  the  equation  for  pressure  thrust  “loss”  gives 
Fpl  =  (Pj  -Pi)  Ab2  =  (6.0969  P,  -  P,)  1 .0972  Abi  =  5.5923  P,  Abi 

which  is  actually  a  thrust  gain.  It  is  a  “gain”  in  thrust,  since  the  direction  of  the  net  force 
due  to  the  pressure  difference  acts  in  the  direction  of  flight  and  acts  to  propel  the  vehicle. 
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Yet  this  does  not  lead  to  something  for  nothing,  since  when  the  momentum  thrust  loss  is 
also  accounted  for,  there  is  a  net  decrease  in  overall  thrust.  This  will  be  shown  in  the 
following  section,  where  the  two  throttling  processes  will  be  compared. 


4.2:  THROTTLING  RESULTS  AND  COMPARISON 

In  the  previous  chapter,  it  was  discovered  that  the  second  solution  to  the  ejector  problem 
pro\  ides  larger  thrust  augmentations,  but  there  were  cases  where  the  entrops’  changes 
were  negative.  So,  to  simplify  presenting  the  trends  of  the  results,  the  second  solution 
will  not  be  covered  in  the  following  comparisons.  Also,  the  entropy  changes  will  not  be 
given,  because  the  entropy  changes  that  occur  from  a  6%  diffuser  loss  do  not  differ 
greatly  from  the  5%  loss  discussed  before. 

After  calculating  the  thrust  decrease  due  to  throttling  the  flow  back  to  ambient  pressure, 
the  following  result  was  achieved. 


Table  8:  Thrust  Augmentation  with  6%  Diffuser 
Loss  and  Bleed  Flow  to  Ambient  Pressure 


Thrust  when  expanded  to  ambient  Pressure  - 
%  Thrust  increase  over  baseline,  -107  kN,  Diffuser  loss=6% 


M//Pr= 

0.90 

0.85 

0.80 

0.75 

0.70 

0.65 

0.60 

0.55 

0.5328 

0.50 

0.45 

0.40 

0.39 

0.49 

0.57 

0.65 

0.73 

0.81 

0.89 

0.97 

1.0000 

1.05 

1.13 

1.20 

0.0 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.1 

0.12 

0.23 

0.25 

0.20 

0.09 

-0.10 

-0.37 

-0.74 

-0.89 

-1.22 

-1.86 

-2.67 

0.2 

0.26 

0.51 

0.62 

0.64 

0.59 

0.46 

0.24 

-0.09 

-0.24 

-0.55 

-1.17 

-1.98 

0.3 

0.36 

0.71 

0.89 

0.96 

0.96 

0.88 

0.69 

0.40 

0.27 

-0.04 

-0.64 

-1.46 

0.4 

0.41 

0.85 

1.07 

1.19 

1.22 

1.17 

1.03 

0.77 

0.64 

0.35 

-0.24 

-1.05 

0.5 

0.42 

0.92 

1.19 

1.33 

1.39 

1.37 

1.26 

1.03 

0.91 

0.64 

0.06 

-0.76 

0.6 

0.39 

0.95 

1.24 

1.41 

1.48 

1.49 

1.41 

1.20 

1.09 

0.83 

0.26 

-0.56 

0.7 

0.32 

0.92 

1.24 

1.42 

1.51 

1.53 

1.47 

1.29 

1.19 

0.95 

0.39 

-0.43 

0.8 

0.21 

0.85 

1.19 

1.39 

1.49 

1.52 

1.48 

1.32 

1.23 

0.99 

0.45 

-0.38 

0.9 

0.07 

0.74 

1.10 

1.31 

1.41 

1.57 

1.46 

1.29 

1.21 

0.98 

0.45 

-0.38 

1.0 

-0.10 

0.60 

0.98 

1.19 

1.29 

1.57 

1.51 

1.21 

1.13 

0.92 

0.39 

-0.44 

1.1 

-0.29 

0.42 

0.81 

1.03 

1.14 

1.52 

1.51 

1.10 

i  1.01 

0.81 

0.30 

-0.54 

1.2 

-0.52 

0.22 

0.62 

0.85 

0.95 

1.43 

1.46 

1.08 

0.86 

0.66 

0.16 

-0.68 

1.3 

-0.76 

-0.01 

0.40 

0.63 

0.74 

1.30 

1.37 

1.02 

0.80 

0.48 

-0.02 

-0.85 

1.4 

-1.03 

-0.26 

0.16 

0.39 

0.50 

1.13 

1.25 

0.92 

0.70 

0.26 

-0.22 

-1.06 

1.5 

-1.32 

-0.54 

-0.11 

0.13 

0.23 

0.93 

1.09 

0.78 

0.57 

0.01 

-0.46 

-1.30 

1.6 

-1.62 

-0.83 

-0.39 

-0.15 

-0.06 

0.70 

0.91 

0.62 

0.41 

-0.15 

-0.72 

-1.56 

1-7 

-1.94 

-1.14 

-0.70 

-0.46 

-0.36 

0.45 

0.70 

0.43 

0.23 

-0.33 

-1.00 

-1.85 
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1,8 

-2.28 

-1.47 

-1.02 

-0.78 

-0.68 

0.17 

0.46 

0.22 

0.01 

-0.54 

1.9 

-2.63 

-1.81 

-1.36 

-1.11 

-1.02 

-0.13 

0.20 

-0.02 

-0.22 

-0.78 

2,0 

-3.00 

-2.17 

-1.71 

-1 .46 

-1.37 

-0.45 

-0.07 

-0.28 

-0.47 

-1.04 

-1.30  -2.15 

-1.63  -2.48 

-1.97  -2.82 


%lncrease  Thrust,  expanded  to  ambient 


Number 


Mass 

Row 

Ratio 


■  1.00-2.00 
□  0.00-1.00 
□ -1.00-0.00 
■ -2.00-1.00 
■ -3.00-2.00 


Figure  16:  Graph  of  %  Thrust  Increase  with 
Bleed  to  Ambient  Pressure 


The  result  of  throttling  the  flow  to  half  of  its  previous  stagnation  pressure  results  in  the 
following  thrust  increases. 


Table  9;  Thrust  Augmentation  with  6%  Diffuser 
Loss  and  Bleed  Flow  to  V2  Po 


Thrust  when  expanded  to  ambient  Pressure  - 
%  Thrust  increase  over  baseline,  -107  kN,  Diffuser  loss=6% 


M^/Pr= 

0.90 

0.85 

0.80 

0.75 

0.70 

0.65 

0.60 

0.55 

0.5328 

0.50 

0.45 

0.40 

0.39 

0.49 

0.57 

0.65 

0.73 

0.81 

0.89 

0.97 

1.0000 

1.05 

1.13 

1.20 

0.0 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.1 

0.64 

0.75 

0.77 

0.72 

0.60 

0.42 

0.15 

-0.22 

-0.38 

-0.71 

-1.34 

-2.15 

0.2 

1.29 

1.54 

1.65 

1.67 

1.62 

1.49 

1.27 

0.94 

0.79 

0.48 

-0.14 

-0.96 

0.3 

1.90 

2.26 

2.43 

2.51 

2.50 

2.42 

2.24 

1.94 

1.81 

1.51 

0.90 

0.09 

0.4 

2.47 

2.91 

3.13 

3.25 

3.28 

3.23 

3.09 

2.82 

2.70 

2.41 

1.82 

1.01 

0.5 

3.00 

3.50 

3.76 

3.90 

3.97 

3.95 

3.84 

3.60 

3.49 

3.21 

2.63 

1.81 

0.6 

3.48 

4.03 

4.33 

4.50 

4.57 

4.58 

4.50 

4.29 

4.18 

3.92 

3.35 

2.53 

0.7 

3.92 

4.52 

4.85 

5.03 

5.12 

5.14 

5.08 

4.90 

4.80 

4.55 

3.99 

3.17 
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0.8 

4.33 

4.97 

5.31 

5.51 

5.60 

5.64 

5.60 

5.44 

5.35 

5.11 

4.57 

3.74 

0.9 

4.71 

5.37 

5.74 

5.94 

6-04 

6.21 

6.09 

5.92 

5.84 

5.62 

5.08 

4.25 

1,0 

5.05 

5.75 

6.12 

6.34 

6.44 

6.72 

6.66 

6.35 

6.28 

6.07 

5.54 

4.71 

1.1 

5.37 

6.09 

6.48 

6.70 

6.80 

7.19 

7.17 

6,76 

6.67 

6.47 

,5.96 

5.13 

1.2 

5.66 

6.40 

6.80 

7.03 

7.13 

7.61 

7.64 

7.26 

7.03 

6.84 

6.34 

5.50 

1.3 

5.93 

6.68 

7.10 

7.33 

7.43 

7.99 

8.06 

7.71 

7.49 

7.17 

6.68 

5.84 

1.4 

6.18 

6.95 

7.37 

7.60 

7.71 

8.34 

8.46 

8.12 

7.91 

7.47 

6.99 

6.15 

1.5 

6.41 

7.19 

7.62 

7.85 

7.95 

8.65 

8.82 

8.51 

8.29 

7.74 

7.27 

6.43 

1.6 

6.62 

7.41 

7.84 

8.08 

8.18 

8.94 

9.15 

8.86 

8.65 

8.09 

7.52 

6.68 

1.7 

6.81 

7.61 

8.05 

8.30 

8.39 

9.20 

9.45 

9.18 

8.98 

8.42 

7.75 

6.91 

1.8 

6.99 

7.80 

8.25 

8.49 

8.58 

9.43 

9.73 

9.48 

9.28 

8.72 

7.96 

7.11 

1.9 

7.15 

7.97 

8.42 

8.67 

8.76 

9.65 

9.99 

9.76 

9.56 

9.00 

8.15 

7.30 

2.0 

7.30 

8.13 

8.59 

8.83 

8.92 

9.84 

10.22 

10.02 

9.82 

9.26 

8.33 

7.47 

%lncrease  Thrust,  expanded  to  ambient 


Bypass 

Mach 

Number 


Mass 

Row 

Ratio 


□  10.00-12.50 

■  7.50-10.00 

□  5.00-7.50 

□  2.50-5.00 

■  0.00-2.50 
■ -2.50-0.00 


Figure  17:  Graph  of  %  Thrust  Increase  with 
Bleed  to  '/2  Po 

Looking  at  the  two  previous  graphs,  it  is  apparent  that  the  process  of  throttling  to  ambient 
pressure  has  a  greater  affect  on  decreasing  realizable  thrust.  When  the  flow  is  only 
throttled  to  V2  Po,  there  is  a  less  drastic  affect  on  overall  thrust.  This  logically  makes 
sense,  since  the  stronger  the  shock  that  the  flow  has  to  pass  through,  the  greater  the 
decrease  in  the  energy  it  has  to  provide  thrust.  The  basic  affect  of  the  throttle  to  ambient 


44 


pressure  is  to  act  like  a  large  drag  brake.  It  almost  completely  removes  any  thrust  gain 
from  the  ejector/engine  combination,  but  it  does  not  completely  destroy  the  thrust  gained. 

Thus,  as  can  be  seen  from  the  first  case,  even  if  the  flow  is  treated  in  a  very  severe 
manner  the  thrust  gained  is  not  completely  destroyed.  Hopefully  this  would  cany  through 
into  any  actual  bleed  configuration  used.  The  only  way  that  the  thrust  could  be 
completely  removed  would  be  to  throttle  to  a  pressure  even  lower  then  ambient,  but  this 
would  be  extremely  irrational.  Hence,  unless  something  illogical  is  done  to  the  flow, 
there  will  still  be  a  thrust  increase  for  a  6%  diffuser  loss.  Larger  diffuser  losses  were  not 
examined,  but  the  thrust  decrease  would  be  taken  in  a  similar  manner  off  of  the  already 
decreased  thrust  profiles  of  the  larger  diffuser  loss  cases. 


CHAPTER  5:  CONCLUSIONS  AND  RECOMMENDATIONS 


5.1:  EJECTORS  AS  A  VIABLE  OPTION? 

In  the  Boeing  design  for  the  HSCT,  ejectors  are  going  to  be  used  for  the  purpose  of  noise 
damping.  In  order  to  determine  if  ejectors  would  provide  thrust  augmentation,  in  the 
supersonic  cruise  configuration,  the  analysis  just  presented  was  conducted. 

5.1.1:  A  Viable  First  Order  Analysis 

The  derivation  of  Chapter  1  (with  results  presented  in  Chapter  2)  showed  a  considerable 
thrust  augmentation.  This  gain  occurred  over  a  wide  range  of  bypass  Mach  numbers  and 
mass  flow  ratios,  when  no  losses,  other  than  mixing  were  taken  into  account. 

Both  the  first  and  the  second  solutions  show  a  large  percentage  increase  over  the  baseline 
thrust.  Looking  at  the  first  solution,  there  are  some  cases  where  the  momentum  equation 
could  not  be  solved,  and  to  reiterate,  these  values  were  boxed  in  to  indicate  that  this  was 
so.  Yet,  even  if  the  boxed  solutions  are  not  possible,  there  is  still  a  large  thrust  increase 
that  can  be  gained. 

As  for  the  second  solution,  the  thrust  increases  are  larger,  but  whether  this  would  be  more 
beneficial  cannot  be  determined  from  the  analysis  presented.  A  model  for  determining 
mixing  length  was  not  included,  so  no  statement  can  be  made  regarding  its  possible 
benefits  over  the  first  solution. 

One  final  thing  to  note  regarding  the  second  solution  is  that  there  are  some  values  of 
b>pass  Mach  number  and  mass  flow  ratio  that  do  not  allow  for  a  solution  because  of  the 
Second  Law  of  Thermodynamics.  The  associated  non-dimensional  entropy  changes  are 
negative,  since  the  second  solution  to  the  momentum  equation  is  supersonic,  but  the 
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incoming  flows  are  both  subsonic.  Therefore,  this  causes  the  negative  entropy  change, 
and  is  not  physically  possible. 

5.1.2:  Effect  of  Diffuser  Loss  on  Thrust 

When  the  affects  of  a  straight  percentage  diffuser  loss  were  applied  to  the  stagnation 
pressure  of  the  secondary  flow,  there  was  a  noticeable  decrease  in  the  thrust 
augmentation.  This  occurred  for  both  the  first  and  the  second  solutions. 

However,  in  both  cases,  there  was  not  a  total  loss  in  thrust  augmentation  until  the  diffuser 
loss  reached  25%.  Therefore,  under  the  assumptions  of  the  analysis  of  this  paper,  ejectors 
still  seem  to  have  promise  even  with  a  diffuser  loss  added. 

5.1.3:  Effect  of  Bleed  Loss  on  Thrust 

Looking  back  at  Chapter  4,  the  bleed  loss  was  modeled  in  two  way.  First,  the  thrust  was 
throttled  directly  to  ambient  pressure,  and  the  second  case  had  the  bleed  flow  lose  half  of 
its  stagnation  pressure.  The  reasons  for  the  second  choice  were  previously,  but  by 
attempting  to  bleed  the  flow  more  “intelligently'’,  the  thrust  losses  experienced  were  not 
nearly  as  large. 

For  the  first  case,  the  losses  were  enough  to  make  any  thrust  gains  from  the  ejector 
practically  disappear.  This  was  mainly  due  to  the  strong  normal  shock  that  was  present 
from  doing  the  bleed  in  the  manner  proposed.  The  throttling  model  that  was  used  created 
such  an  adverse  affect  on  mass  flow;  the  result  was  an  essential  drag  brake. 

But.  doing  the  bleed  so  that  there  was  only  a  loss  in  half  the  stagnation  pressure,  the 
results  are  more  tolerable,  and  all  of  the  thrust  augmentation  is  not  destroyed. 
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5.2:  WHAT  NEEDS  TO  BE  DONE 

After  doing  this  analysis,  the  author  determined  ejectors  seem  very’  promising  for  further 
research.  Some  of  the  areas  that  can  be  looked  into  are.  extending  the  mathematical 
modeling  of  the  problem  itself,  and  beginning  experimental  testing. 

5.2.1:  The  Mathematical  Model 

Some  of  the  shortcomings  of  the  model  used  in  the  analysis  were,  for  example;  all  gases 
were  assumed  to  be  ideal,  and  the  ratio  of  specific  heats  was  constant.  While  these  make 
good  first  order  assumptions,  the  favorable  results  achieved  thus  far  lead  the  author  to 
conclude  that  the  model  should  be  refined  and  examined  further. 

One  of  the  first  adaptations  to  make  is  to  longer  assume  that  the  gases  that  will  be  mixing 
are  ideal.  Real  gas  effects  can  be  added  into  the  model,  and  their  effect  on  the  analysis 
may  then  be  determined. 

Second,  the  gases  were  all  assumed  to  be  air.  and  this  would  not  be  so  in  the  actual 
working  model  for  the  HSCT.  There  would  be  combustible  gases  expelled  from  the 
primary  core  (the  engine)  and  these  could  possibly  combust  with  the  secondary  (entrained 
air)  flow.  So,  adding  combustion  to  the  model  seems  a  reasonable  extension. 

Third,  the  author's  analysis  was  essentially  an  analysis  of  the  performance  of  an  engine 
that  was  uninstalled.  No  external  drag  forces  were  accounted  for,  so  if  desired,  one  could 
examine  in  greater  detail  the  effect  of  drag  on  the  system  over  various  possible  inlet 
geometry’s. 

Fourth,  there  is  no  model  being  used  to  determine  the  exact  method  of  mixing  in  the 
thrust  (mixing)  chamber.  It  was  assumed  that  mixing  was  complete  when  the  mixing 
Mach  number  (A/^)  was  determined.  But.  nothing  was  done  to  calculate  the  distance  that 
would  be  required  to  mix.  So,  a  very  interesting  area  of  research  would  be  an 
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examination  of  how  to  make  the  two  flows  mix  in  as  short  a  distance  as  possible.  This 
would  make  the  thrust  augmentation  an  even  greater  benefit,  since  the  weight  of  the 
ejector  assembly  used  to  generate  thrust  would  decrease  simply  due  to  the  shortened 
mixing  length. 

Finally,  a  different  model  for  diffuser  and  bleed  loss  could  be  attempted.  The  exact 
manner  in  which  it  was  done  by  the  author,  may  not  be  the  most  elegant  way  to  do  so. 
Solving  the  exact  characteristics  of  the  inlet  diffuser  loss,  and  bleed  loss,  and  how  to 
minimize  both  is  an  area  of  significant  importance  to  the  success  of  this  endeavor.  If 
these  two  could  be  kept  as  low  as  possible,  the  thrust  augmentation  achieved  is  of  e\  en 
greater  value. 

5.2.2;  Experimental  Testing 

Above  and  beyond  the  intricacies  of  theoretical  modeling,  some  actual  hardware  analyses 
and  experiments  should  be  conducted.  Empirical  testing  sometimes  serves  as  a  guide  to 
the  theoretician,  and  would  be  very  helpful  for  this  project.  The  exact  dynamic  that 
causes  the  solution  of  mixing  to  literally  choke,  and  sometimes  be  impossible  is  not 
readily  apparent  from  the  math  of  theory  alone.  Examining  this  phenomena,  and  the 
thrust  augmentations  that  appear  from  theory,  would  provide  valuable  data  for  production 
of  the  propulsion  system/ejector  combination  for  the  HSCT. 


5.3:  FINAL  THOUGHTS 

Looking  back  on  everything  that  has  been  done  so  far,  ejectors  seem  very  promising  for 
use  as  a  thrust  augmentation  device  on  the  HSCT  during  supersonic  cruise.  Because  this 
was  only  a  first  order  analysis,  more  can  be  done,  and  more  should  be  done!  Ejectors  are 
already  proposed  for  use  during  the  takeoff  phase,  since  they  can  be  used  for  noise 
damping,  so  their  possible  use  in  the  cruise  configuration  should  not  be  ignored. 
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It  is  the  author's  recommendation  that  this  problem  be  studied  further  in  some  if  not  all  of 
the  ways  mentioned  above.  Then  the  viability  of  using  ejectors  as  a  thrust  augmenting 
device  in  supersonic  cruise  can  be  determined  definitively  one  way  or  the  other. 


ENDNOTES 


1 .  See  Appendix  A  for  the  data  provided  by  the  Boeing  Company. 

2.  See  page  563  of  Introduction  to  Flight.  3^*^.  by  Anderson  for  standard  atmosphere  data. 

3.  See  pages  134-5  of  Liepmann  and  Roshko  for  continuity,  momentum  and  energy 
equation  derivations. 

4.  See  Appendix  B  for  an  analytic  solution  to  Me  as  derived  by  Alperin  and  Wu  from 
Thrust  Augmenting  Ejectors,  Part  I. 

5.  See  page  1 89-192  of  Aerothermodynamics  ...  by  Oates  for  the  derivation  of  the  thrust 
equation. 

6.  See  page  123  of  Physical  Gas  Dynamics,  by  Vincenti  and  Krueger. 

7.  See  page  33  of  Gas  Dynamics  Tables  for  Air  by  Emmons. 
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APPENDIX  A:  BOEING  ENGINE  DATA 


According  to  the  Boeing  engine  data  given,  the  relevant  information  for 
the  primar\’  stream  and  the  ambient  conditions  are,  altitude,  flight  mach  number,  ambient 
pressure,  ambient  temperature,  true  airspeed,  core  stream  total  temperature,  core  stream 
total  gas  flow  at  mixing  plane,  and  core  stream  total  pressure  at  mixing  plane. 

The  conditions  given  were: 


Alt  =  50.000  ft 
Mach  Number  =  2.4 
Ambient  Pressure  =  1 .682  psia 
Ambient  Temperature  =  390  R 
True  Airspeed  =  2324  ft/s 
Total  Temperature  =  2277  R 
Total  Gas  Flow  =  436.92  Ibm/s 
Total  Pressure  =  46.96  psia 

All  of  the  data,  except  for  altitude,  will  be  converted  to  SI  units  for  ease  of 
calculation.  Since  altitude  is  not  a  calculation  variable,  it  has  not  been  converted,  and  it 
also  gives  a  better  intuitive  feel  as  to  the  flight  level  in  most  aviation  related  fields. 
Knowing  that  1  atm  =  14.7  psia  =  101,350  Pa  allows  for  easy  conversion  of  the  pressure 
terms,  and  temperature  is  easily  converted  since  1°  K  =  1.8°  R.  Other  useful  conversion 
factors  are  1  ft  =  0.3048  m  and  1  Ibm  =  0.454  kg.  Thus  converting  these  values  gives: 

Alt  =  50,000  ft 
Mach  Number  =  2.4 
Ambient  Pressure  =  1 1 ,597  Pa 
Ambient  Temperature  =  216.7  K 
True  Airspeed  «  708  m/s 
Total  Temperature  =  1265  K 
Total  Pressure  =  323,780  Pa 

While  this  data  comes  from  the  Boeing  Company,  the  values  of  ambient 
pressure  and  temperature  were  cross  checked  with  tables  of  a  standard  atmosphere  and 
the  value  for  true  airspeed  was  calculated  directly  from  the  definition  of  Mach  number 
and  the  given  value  of  temperature.  Also,  the  values  of  total  temperature  and  pressure 
that  were  given  for  the  core  flow  had  to  be  used  to  model  the  engine  they  used. 

Here  is  the  data  in  raw  form: 
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Table  10:  Boeing  Data  Format 


ALT 

MN 

DTAMB 

XXX 

XXX 

PAMR 

TAMB 

m 

XXX 

XXX 

XXX 

XXX 

XXX 

XXX 

XXX 

XXX 

XXX 

TT6 

XXX 

XXX 

XXX 

XXX 

XXX 

XXX 

XXX 

WG8 

TTS 

PTS 

XXX 

XXX 

ASGEOSTW 

A9STW 

WGMIX 

TTMIX 

PTMIX 

XXX 

XXX 

A9MIX 

XXX 

VJI 

XXX 

XXX 

XXX 

XXX 

XXX 

XXX 

XXX 

XXX 

XXX 

XXX 

XXX 

XXX 

XXX 

XXX 

XXX 

XXX 

XXX 

XXX 

XXX 

XXX 

NPR 

ASGEOSP 

A9SF 

XXX 

TT16 

PT16 

WG16 

WG5  5 

XXX 

XXX 

XXX 

XXX 

XXX 

XXX 

T155. 

TT7 

XXX 

XXX 

XXX 

XXX 

XXX 

XXX 

PTRR 

XXX 

WS 

PTS 

TTS 

VJSBP 

AJSESP 

WSMIX 

PTMIX 

TTMIX 

VJIM 

A9MIXSP 

Table  1 1 :  Boeing  Raw  Data 

50000  ■ 

2 .40 

.0 

00 . 00 

.  0 

i.6a2 

390.0 

232.^ 

0 . 00 

00 . 00 

000 . 0 

000 . 00 

000 . 00 

.  0 

000 . 0 

0000  . 

000 . 0 

2277  . 

00000  . 

00000  . 

0 .000 

.  0 

00 . 00 

0000  . 

0000 . 0 

701 . 94 

1S52  . 

46.45 

0 . 000 

000 . 0 

1290  . 

795  . 

.  0 

.  0 

.0 

.  0 

.  0 

.  0 

00000 . 

754  . 

0 . 000 

.0 

.0 

00000  . 

000 . 00 

0 . 000 

0 . 000 

.  00 

0 . 00 

0 . 00 

0 . 000 

0000. 

000. 

0 . 000 

00 . 00 

0 . 000 

.000 

0000 

00000,. 

000 . 00 

27.62 

8 . 958 

3 .30 

00 . 00 

1079  . 

49.31 

265 . 02 

436 ■ 92 

.  000 

0 . 000 

0.000 

0000  . 

0000  . 

000 . 0 

2277. 

1852  . 

000  . 

0000 . 0 

000 . 0 

0 . 0000 

0 .000 

000 . 0 

4£...a£ 

0 . 000 

The  relevant  data  in  the  above  tables  are  underlined,  and  were  converted  on  the  previot 

page.  The  underlined  abbreviations  mean  in  the  previous  table  mean: 


ALT  -  Geopotential  Altitude  -  FT 
MN  -  Flight  Mach  Number 
PAMB  -  Ambient  Pressure  -  PSIA 
TAMB  -  Ambient  Temperature  -  R 
VO  -  True  Airspeed  -  FT/SEC 

WG55  -  Core  Stream  Total  Gas  Flow  at  Mixing  Plane  -  LBM/SEC 
TT55  -  Core  Stream  Total  Temperature  at  Mixing  Plane  -  R 
PT55  -  Core  Stream  Total  Pressure  at  Mixing  Plane  -  PSIA 


APPENDIX  B:  ALPERINAVU  SOLUTION  TO  THE  EJECTOR  PROBLEM 


In  the  original  AlperinAVu  derivation  in  “High  Speed  Ejectors”  or  in  “Thrust 
Augmenting  Ejectors,  Part  I”,  different  terminology  was  used  for  their  station  locations. 
In  this  appendix,  I  will  do  the  derivation  in  the  terms  of  the  stations  previously  defined  in 
the  main  body  of  the  text. 

For  Alperin  and  Wu,  they  used  the  subscripts  i,  p  and  2  to  respectively  describe 
the  ejector,  the  core  and  the  mixing  plane.  For  the  purpose  of  this  appendix,  the  station 
numbers  5,  4  and  6  will  be  used  in  place  of  the  subscripts  just  described.  So,  just  to  be 
clear,  5  =  ejector,  4  =  primary  (core)  and  6  =  mixed  flow. 

Alperin  and  Wu’s  derivation  clearly  shows  the  presence  of  the  two  solutions  to 
the  mixing  plane  Mach  number.  It  also  gives  criteria  for  when  it  is  possible  that  there 
will  only  be  one  root.  For  their  derivation  the  following  assumptions  were  made: 

1)  All  fluids  are  compressible  and  thermally  and  calorically  perfect. 

2)  Skin  friction  and  blockage  losses  are  neglected. 

3)  Mixing  is  initiated  in  a  constant  area  duct  at  the  location  where  the  primary  flow  is 
fully  expanded  (primary  flow  pressure  is  equal  to  the  local  secondary  flow  pressure  P4  = 
P5). 

4)  Complete  mixing  occurs  in  a  constant  cross-sectional  channel. 

5)  Ejector  surfaces  are  adiabatic. 

Clearly,  these  are  very  similar  to  the  assumptions  previously  made  in  the  main 
body,  and  so  should  lead  to  similar  results.  To  begin,  several  beisic  equations  need  to  be 
defined.  They  are: 

Continuity: 

•  •  • 

mb  =  ms  +  m\  —  p  i  Ai  V5  +  p  A  Aa  VA  =  p  b  Ab  V6  (Bl) 

Energy: 


mA  Toa  +  ms  Tos  —  (mA  +  ms)  T06 


(B2) 
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Momentum: 

•  •  •  • 

(Ps  -  P6)A6  +  ms  V5  +OT4  -  {ms  +tm)V6  =  0 
Area  Ratios: 
a  =  A(,IA4 

and  since  As  =  Aj+A4 

(a-l)  =  A,/A4 

Stagnation  Temperature  Relations: 

Tos  =  75(^1  + 

Toe  =  76(^1  + 


7'o4  =  n  (^1  +  Ma- 


Acoustic  Velocity: 

V5  =  Msyjy  R  Ts 

V6  =  RTb 

V4  =  MA-yj/  RTa 

Static  Pressure: 

mA  R  Ta 


(B3) 

(B4) 

(B5) 

(B6) 

(B7) 

(B8) 

(B9) 

(BIO) 

(Bll) 


Ps  = 


VA  As 


=  Pa 


(B12) 
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p  _  {ms  +mi)RTb 
V6  A  6 

Static  Temperature: 

Tos 

r-  1 
2 


Tb  = 


Tob 


Ta 


7’o4 


(B13) 


(B14) 


(B15) 


(B16) 


Mass  Flow  Ratios: 


It  is  easy  to  prove  the  right  hand  side  of  equation  (B17),  if  one  takes  the  first 
equation  after  the  variable  and  substitutes  in  the  definition  of  mass  flow.  This  gives 

_  psAsvs  _  RTaPs 
^  p  A  Aa  V A  RTs  Pa 

The  right  hand  side  comes  from  the  application  of  the  perfect  gas  law,  (B5),  (B9) 
and  (Bll).  Canceling  similar  variables  from  top  and  bottom,  remembering  the  primary 
flow  is  fully  expanded  so  P\  =  P$,  and  applying  equations  (B16)  and  (B17)  one  has 
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Now,  with  the  basic  definitions,  one  can  take  the  momentum  equation  above  (B3) 
and  by  expressing  pressures  in  terms  of  the  perfect  gas  law,  velocity  in  terms  of  Mach 
number  and  temperature  in  terms  of  stagnation  temperature,  iV/is  can  be  solved  for  directly. 

Taking  (B3)  and  substituting  (B12)  and  (B13)  gives: 

•  •  •  • 

Ae  +  ms  vs  +OT4  im  -  (ms  +/W4)v6  =  0 


Dividing  through  by  m4  and  multiplying  the  back  into  the  parentheses,  gives  (in  view 
of  equations  (B1 7)  and  (B4)) 


RT4a  (I  +  /j)RT6  ^  ^  ^  .  -A 

-  -  -  +  yU  V5  +  V4  -  (1  +  /i)  V6  —  0 


Substituting  (B9)  -  (B1 1)  into  the  above,  then  canceling  the  appropriate  values  of  T  and  R 
gives 

a  yjR  Tj  _  (1  +  fx)^lRT6  fxMsJ/RTs  +  M4^yRT4  -  (1  +//)  Mb^/RTh  =  0 

{M4^  A/5  77  j 

Multiplying  by  l/ yj R  y  and  applying  equations  (B14)  -  (B16)  gives 
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Gathering  like  terms  and  multiplying  by 


-(1  +  //)'v/7'06 


gives 


1+/  Me- -■ 


ay  Me  1 

iToi 

iToi 

(l  +  ju)  y  Mi  V 

706  1 

iV'.v/r  (>^>')  ^ 

'ro6i 

1  +  ^  Me- 


1  +  ^M4- 


//  /  Ms  Me  To5  To4 

(i+/y)  vr^v^ 


1  +  ^  Me- 


1  + 


/-I 


=  0 


Ms- 


The  last  three  terms  on  the  left  hand  side  all  have 
factoring  it  out  gives 


-r 

(1  +  //) 


r-i 

7 


Me- 


1  +  /  Me' 


(1  +  /^) 


7*06 


ITos 

/uMi 

Mi 

i  foe 

1 

vro4 

1 

\  +  Mi- 

V  2 

1+- 


a 


y  ) 


y  X/6^1  +  M6“  =  0  (equation  continued  from  previous  line) 

so  letting 


J  = 


1  Toa 


(i+/i)  vro6 


ITos 

/u  Ms 

J  Toe 

7o5  Ma 


-f - 


1  + 


a 


one  has 


\+y  Me-  -Jy  Me^l  +  ^  Me'  =  0 


,  SO 
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moving  the  third  term  of  the  left  hand  side  to  the  right  hand  side  and  squaring  both  sides 
gives 


1+2/  Mb-  +(/  Me-)'  =J-y- Mb' 


1  +  ^  Mb- 
2 


multiplying  through  on  the  right  hand  side  and  letting  X=  y  M^  gives 
1  +  2X  +X-  =  yf-X  +  X-J-  ^ 


gathering  like  terms  in  X  gives 
1+  {2-y  J-]X  +{\-J-  =  0 


(B18) 


then  using  the  quadratic  formula  to  solve  forX  one  gets 

^  -B  ±  Vb-  -4  a  ^  ,2  j  A  ,  r2  /-J 

X  =  - if  B  =  2-y  J  and  A  =  1-  J  - 

2  A  2 

knowing  X=  y  M^,  solving  for  Mt  from  the  above  gives 


Mb 


-B  ±Vb-  -4  a 


2  y  A 


(B19) 


However,  the  above  derivation  leaves  on  thing  to  question,  what  is  Tog? 
According  to  Alperin  and  Wu 


Tob 

TZ 


1  +  fi 


To4 


1+// 


(B20) 


This  is  easily  derived  from  (B2).  Dividing  the  left  and  right  by  To-t  {mA+  ms) ;  then 
moving  switching  the  order  of  the  sides  gives 


1 


Toe 


(/774+  mi) 


7W4  +  ms 


tJ 


then  factoring  out  primary  mass  flow  {mi)  and  rewriting  everything  in  terms  of  mass 
flow  ratio  (/^ )  gives  (B20). 

So,  for  any  given  set  of  flow  properties  {Ms,  M4,  T04,  T05)  and  area  ratio  {a),  there 
are  two  possible  solutions  to  (B19).  Alperin  and  Wu  refer  to  these  Mach  numbers  as  Msn) 
(Me.first)  when  the  negative  sign  of  (B19)  is  used  and  .V/2(ii)  (Me.second)  when  the  positive 
sign  is  used.  These  are  the  first  and  second  solutions  that  were  referred  to  in  the  main 
paper. 


As  Alperin  and  Wu  note  in  their  paper,  it  is  obvious  that  Me  has  no  real  value  if  B‘ 
<  4A  and  has  two  real  solutions  if  B^  >  4A.  These  two  real  solutions  are  the  subsonic  and 
supersonic  solutions  to  the  momentum  equation  in  the  mixing  plane.  Alperin  and  Wu 
state  the  following  relationship  between  the  two  Mach  numbers  in  “Thrust  Augmenting 
Ejectors,  Part  I”,  but  they  prove  a  similar  result  in  and  appendix  of  “High  Speed 
Ejectors”. 


= 


+  2 

2}'  M2{\\)^  ■  0 


(B21) 


This  is  the  relationship  between  Mach  numbers  across  a  normal  shock  wave. 
However,  when  B^  =  4A,  Me  has  only  one  real  value,  and  at  this  condition,  A/e.nrst  = 

jV/ij  second  ~  1 »  sinCC 

J  =  Jt  =  (1  /  y).yj2[}'  +  l)  when  B^  =  4A  (B22) 

This  is  the  mechanical  choking  which  was  mentioned  in  the  main  paper;  however, 
it  is  not  directly  seen  in  any  of  the  cases  presented.  Its  occurrence  in  terms  of  the  main 
section  of  the  paper  would  be  indicated  by  a  end  of  mixing  Mach  number  equal  to  1  and  a 
solver  solution  equal  to  zero.  Therefore  if  J  =  Jc,  when  B^  =  4A,  the  flow  will  choke  and 
there  will  be  only  one  solution. 

The  combination  of  factors  to  give  B^  =  4A  was  never  shown  in  the  cases  ran  in 
the  creation  of  the  results  in  the  main  paper,  but  it  could  realistically  occur.  This  would 
essentially  be  the  case  if  the  supersonic  flow  of  the  primary  stream  expanded  to  Mach  1 
while  the  secondary  stream  was  “virtually”  compressed  by  the  primary  flow  and  therefore 
accelerated  from  a  subsonic  value  up  to  Mach  1 .  This  would  result  in  a  choked  solution. 
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One  final  comment  regarding  the  Alperin  and  Wu  solution  as  compared  to  the 
method  used  in  the  main  section  is  that  they  chose  to  fix  A  and  then  var\'  different 
parameters,  yet,  the  results  are  essentially  the  same  as  the  results  achieved  by  the  analysis 
conducted  by  the  main  paper.  The  results  are  not  directly  comparable,  but  a  point  by 
point  comparison  can  be  made,  and  the  two  results  are  found  to  match. 


APPENDIX  C:  DIFFUSER  LOSS  DATA 


5%  Diffuser  Loss  Data: 


First  Solution; 


Table  12:  Table  of  Entropy  Change  5%  Diffuser 
Loss  ( 1  St  Solution) 


Values  of 

ASinl 

Pr  =  0.9003 

Pr  =  0.70 

Pr  =  0.55 

Pr  =  0.4 

• 

mi  R 

//=0.0 

2.36E‘°’ 

-5.6E-'^ 

-l.OE''-^ 

-3.8E-'’’ 

//=  1.0 

0.939361 

0.872174 

0.876349 

0.957274 

//=2.0 

1.37646 

1.294843 

1.23408 

1.369679 

Table  13:  Thrust  Augmentation  with  5%  Diffuser 
Loss  (1st  Solution) 


Thrust  when  expanded  to  ambient  Pressure  - 
%  Thrust  increase  over  baseline,  -107  kN 


M//Pr= 

0.90 

0.85 

0.80 

0.75 

0.70 

0.65 

0.60 

0.55 

0.5328 

0.50 

0.45 

0.40 

0.39 

0.49 

0.57 

0.65 

0.73 

0.81 

0.89 

0.97 

1.0000 

1.05 

1.13 

1.22 

0.0 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.1 

0.68 

0.79 

0.81 

0.77 

0.66 

0.48 

0.21 

-0.15 

-0.30 

-0.63 

-1.25 

-2.06 

0.2 

1.38 

1.62 

1.72 

1.74 

1.70 

1.57 

1.35 

1.03 

0.89 

0.58 

-0.03 

-0.84 

0.3 

2.02 

2.36 

2.53 

2.60 

2.60 

2.52 

2.35 

2.06 

1.93 

1.63 

1.03 

0.22 

0.4 

2.61 

3.04 

3.26 

3.37 

3.40 

3.36 

3.22 

2.96 

2.84 

2.55 

1.97 

1.16 

0.5 

3.17 

3.65 

3.91 

4.05 

4.11 

4.09 

3.98 

3.75 

3.64 

3.37 

2.80 

1.99 

0.6 

3.67 

4.22 

4.50 

4.66 

4.74 

4.74 

4.66 

4.46 

4.35 

4.10 

3.54 

2.73 

0.7 

4.14 

4.73 

5.04 

5.21 

5.30 

5.32 

5.26 

5.09 

4.99 

4.75 

4.20 

3.38 

0.8 

4.58 

5.20 

5.53 

5.71 

5.81 

5.89 

5.80 

5.65 

5.56 

5.33 

4.79 

3.97 

0.9 

4.97 

5.62 

5.98 

6.17 

6.26 

6.48 

6.37 

6.15 

6.07 

5.85 

5.32 

4.50 

1.0 

5.34 

6.02 

6.38 

6.59 

6.68 

7.02 

6.95 

6.60 

6.53 

6.32 

5.80 

4.98 

1.1 

5.68 

6.38 

6.76 

6.97 

7.09 

7.50 

7.49 

7.09 

6.94 

6.74 

6.24 

5.42 

1.2 

6.00 

6.71 

7.10 

7.32 

7.46 

7.94 

7.97 

7.60 

7.38 

7.13 

6.63 

5.81 

1.3 

6.29 

7.02 

7.42 

7.64 

7.79 

8.35 

8.42 

8.07 

7.85 

7.48 

6.99 

6.17 

1.4 

6.56 

7.30 

7.71 

7.93 

8.08 

8.71 

8.83 

8.50 

8.29 

7.79 

7.32 

6.49 

1.5 

6.81 

7.57 

7.98 

8.21 

8.35 

9.05 

9.21 

8.91 

8.69 

8.14 

7.62 

6.79 

1.6 

7.04 

7.81 

8.23 

8.46 

8.58 

9.35 

9.56 

9.28 

9.07 

8.52 

7.90 

7.06 

1.7 

7.25 

8.03 

8.46 

8.69 

8.80 

9.63 

9.89 

9.62 

9.42 

8-86 

8.15 

7.31 

1.8 

7.45 

8.24 

8.67 

8.90 

8.98 

9.89 

10.18 

9.94 

9.74 

9.18 

8.37 

7.54 

1.9 

7.63 

8.43 

8.87 

9.10 

9.18 

10.12 

10.46 

10.24 

10.04 

9.48 

8.58 

7.74 

2.0 

7.80 

8.61 

9.05 

9.29 

9.36 

10.33 

10.72 

10.52 

10.32 

9.76 

8.78 

7.93 
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Second  Solution: 


Table  14:  Table  of  Entropy  Change  5%  Diffuser 
Loss  (2nd  Solution) 


Values  of 

AiSrof 

Pr  =  0.9003 

Pr  =  0.70 

Pr  =  0.55 

II 

o 

• 

m^  R 

o 

o 

II 

^  tr*! 6 

-j.jh 

-5.6E‘'‘^ 

-l.OE''-' 

-5.6E-'’ 

1.0 

0.303354 

0.87217 

0.876095 

0.906265 

o 

CN 

II 

-0.64911 

1.294605 

1.23408 

1.344005 

Table  15:  Thrust  Augmentation  with  5%  Diffuser 
Loss  (2nd  Solution) 


Thrust  when  expanded  to  ambient  Pressure  - 
%  Thrust  increase  over  baseline,  -107  kN 


M^/Pr= 

0.90 

0.85 

0.80 

0.75 

0.70 

0.65 

0.60 

0.55 

0.5328 

0.50 

0.45 

0.40 

0.39 

0.49 

0.57 

0.65 

0.73 

0.81 

0.89 

0.97 

1.0000 

1.05 

1.13 

1.22 

0.0 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.1 

1.13 

1.04 

1.03 

1.04 

1.04 

1.03 

1.02 

1.01 

1.00 

0.98 

0.95 

0.91 

0.2 

2.46 

2.05 

1.97 

1.96 

1.96 

1.95 

1.94 

1.91 

1.89 

1.86 

1.80 

1.73 

0.3 

3.95 

3.03 

2.82 

2.78 

2.78 

2.77 

2.75 

2.71 

2.69 

2.65 

2.57 

2.46 

0.4 

5.56 

4.00 

3.61 

3.52 

3.52 

3.51 

3.49 

3.44 

3.41 

3.36 

3.26 

3.13 

0.5 

7.29 

4.96 

4.34 

4.19 

4.18 

4.17 

4.15 

4.09 

4.07 

4.01 

3.89 

3.73 

0.6 

9.13 

5.91 

5.03 

4.81 

4.78 

4.78 

4.75 

4.69 

4.66 

4.59 

4.45 

4.27 

0.7 

11.06 

6.86 

5.68 

5.36 

5.32 

5.33 

5.30 

5.23 

5.20 

5.12 

4.97 

4.77 

0.8 

13.07 

7.80 

6.29 

5.88 

5.82 

5.89 

5.80 

5.73 

5.69 

5.61 

5.44 

5.22 

0.9 

15.16 

8.75 

6.88 

6.35 

6.27 

6.48 

6.37 

6.18 

6.14 

6.05 

5.87 

5.63 

1.0 

17.33 

9.70 

7.45 

6.79 

6.68 

7.02 

6.95 

6.60 

6.56 

6.46 

6.26 

6.00 

1.1 

19.56 

10.65 

8.00 

7.20 

7.09 

7.50 

7.49 

7.09 

6.94 

6.83 

6.62 

6.35 

1.2 

21.86 

11.60 

8.53 

7.59 

7.46 

7.94 

7.97 

7.60 

7.38 

7.18 

6.95 

6.66 

1.3 

24.22 

12.56 

9.05 

7.95 

7.79 

8.35 

8.42 

8.07 

7.85 

7.50 

7.26 

6.95 

1.4 

26.63 

13.53 

9.55 

8.29 

8.08 

8.71 

8.83 

8.50 

8.29 

7.79 

7.54 

7.22 

1.5 

29.09 

14.50 

10.05 

8.61 

8.35 

9.05 

9.21 

8.91 

8.69 

8.14 

7.80 

7.46 

1.6 

31.60 

15.48 

10.53 

8.92 

8.58 

9.35 

9.56 

9.28 

9.07 

8.52 

8.04 

7.69 

1.7 

34.16 

16.47 

11.01 

9.21 

8.80 

9.63 

9.89 

9.62 

9.42 

8.86 

8.26 

7.90 

1.8 

36.76 

17.46 

11.48 

9.49 

8.98 

9.89 

10.18 

9.94 

9.74 

9.18 

8.46 

8.09 

1.9 

39.40 

18.45 

11.95 

9.76 

9.18 

10.12 

10.46 

10.24 

10.04 

9.48 

8.65 

8.26 

2.0 

42.07 

19.46 

12.41 

10.02 

9.37 

10.33 

10.72 

10.52 

10.32 

9.76 

8.83 

8.43 
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10%  Diffuser  Loss  Data: 


First  Solution: 


Table  16:  Table  of  Entropy  Change  10%  Diffuser 
Loss  (1st  Solution) 


Values  of 
ASioi 

Pr  =  0.9003 

Pr  =  0.70 

Pr  =  0.55 

II 

o 

• 

mi  R 

o 

o 

II 

7.8  lE’^^. 

5.55E-'^ 

-4.3E‘" 

JU=  1.0 

0.958441 

0.879811 

0.884675 

0.97076 

II 

b 

1.398703 

1.302486 

1.257077 

1.380831 

Table  17:  Thrust  Augmentation  with  10% 
Diffuser  Loss  (1st  Solution) 


Thrust  when  expanded  to  ambient  Pressure 
%  Thrust  increase  over  baseline,  -107  kN 


^!Pr= 

0.90 

0.85 

0.80 

0.75 

0.70 

0.65 

0.60 

0.55 

0.5328 

0.50 

0.45 

0.40 

0.39 

0.49 

0.57 

0.65 

0.73 

0.81 

0.89 

0.97 

1.0000 

1.05 

1.13 

1.22 

0.0 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.1 

0.43 

0.54 

0.56 

0.50 

0.37 

0.16 

-0.14 

-0.53 

-0.69 

-1.05 

-1.71 

-2.55 

0.2 

0.94 

1.22 

1.33 

1.36 

1.29 

1.15 

0.90 

0.55 

0-39 

0.06 

-0.59 

-1.44 

0.3 

1.43 

1.82 

2.02 

2.10 

2.09 

1.99 

1.79 

1.47 

1.33 

1.01 

0.38 

-0.47 

0.4 

1.88 

2.37 

2.62 

2.75 

2.78 

2.73 

2.56 

2.28 

2.15 

1.84 

1.22 

0.37 

0.5 

2.30 

2.86 

3.15 

3.32 

3.38 

3.36 

3.23 

2.98 

2.85 

2.56 

1.95 

1.10 

0.6 

2.67 

3.29 

3.63 

3.82 

3.91 

3.91 

3.82 

3.59 

3.47 

3.19 

2.59 

1.74 

0.7 

3.02 

3.69 

4.05 

4.26 

4.37 

4.40 

4.32 

4.12 

4.01 

3.75 

3.16 

2.30 

0.8 

3.33 

4.04 

4.43 

4.66 

4.78 

4.82 

4.77 

4.59 

4.49 

4.23 

3.66 

2.80 

0.9 

3.61 

4.35 

4.77 

5.01 

5.14 

5.18 

5.15 

5.00 

4.90 

4.66 

4.09 

3.23 

1.0 

3.87 

4.64 

5.07 

5.32 

5.45 

5.51 

5.49 

5.36 

5.27 

5.04 

4.48 

3.62 

1.1 

4.09 

4.89 

5.33 

5.59 

5.73 

5.90 

5.87 

5.67 

5.59 

5.37 

4.82 

3.95 

1.2 

4.30 

5.11 

5-57 

5.84 

5.98 

6.24 

6.26 

5.94 

5.87 

5.66 

5.12 

4.25 

1.3 

4.48 

5.32 

5.79 

6.06 

6.20 

6.54 

6.61 

6.24 

6.12 

5.92 

5.39 

4.52 

1.4 

4.64 

5.49 

5.97 

6.25 

6.39 

6.81 

6.93 

6.58 

6.36 

6.14 

5.63 

4.75 

1.5 

4.79 

5.65 

6.14 

6-42 

6.57 

7.05 

7.21 

6-89 

6.67 

6.34 

5.83 

4.95 

1.6 

4.91 

5.79 

6.29 

S.57\2 

6.48 

7.25 

7.46 

7.16 

6.95 

6.51 

6.01 

5.13 

1.7 

5.02 

5.91 

6.42 

6.70 

6.85 

7.44 

7.69 

7.41 

7.20 

6.66 

6.17 

5.28 

1.8 

5.12 

6.02 

6.53 

6.82 

6.96 

7.59 

7.89 

7.64 

7.43 

6.85 

6.31 

5.41 

1.9 

5.20 

6.11 

6.62 

6.92 

7.06 

7.73 

8.07 

7.84 

7.63 

7.06 

6.42 

5.53 

2.0 

5.27 

6.19 

6.71 

7.00 

7.14 

7.85 

8.23 

8.02 

7.82 

7.24 

6.52 

5.62 
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Second  Solution: 


Table  18:  Table  of  Entropy  Change  10%  Diffuser 
Loss  (2nd  Solution) 


Values  of 

ASini 

Pr  =  0.9003 

Pr  =  0.70 

Pr  =  0.55 

Pr  =  0.4 

• 

7774  R 

//=0.0 

2.78E-'^ 

5.55E-'^ 

LllE*'^ 

//=  1.0 

0.303683 

0.877634 

0.881424 

0.909129 

//=2.0 

-0.65814 

1.29975 

1.257077 

1.347212 

Table  19:  Thrust  Augmentation  with  10% 
Diffuser  Loss  (2nd  Solution) 


Thrust  when  expanded  to  ambient  Pressure  - 
%  Thrust  increase  over  baseline,  -107  kN 


M//Pr= 

0.90 

0.85 

0.80 

0.75 

0.70 

0.65 

0.60 

0.55 

0.5328 

0.50 

0.45 

0.40 

0.39 

0.49 

0.57 

0.65 

0.73 

0.81 

0.89 

0.97 

1.0000 

1.05 

1.13 

1.22 

0.0 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.1 

0.98 

0.90 

0.89 

0.90 

0.90 

0.90 

0.89 

0.88 

0.87 

0.86 

0.82 

0.79 

0.2 

2.18 

1.78 

1.70 

1.69 

1.69 

1.69 

1.68 

1.65 

1.64 

1.61 

1.55 

1.48 

0.3 

3.55 

2.65 

2.43 

2.39 

2.39 

2.39 

2.37 

2.33 

2.32 

2.28 

2.20 

2.10 

0.4 

5.07 

3.50 

3.10 

3.01 

3.00 

3-00 

2.98 

2.94 

2.92 

2.87 

2.77 

2.64 

0.5 

6.70 

4.35 

3.72 

3.57 

3.55 

3.55 

3.53 

3.48 

3.45 

3.40 

3.28 

3.13 

0.6 

8.45 

5.20 

4.30 

4.07 

4.03 

4.03 

4.01 

3.96 

3.93 

3.87 

3.74 

3.56 

0.7 

10.30 

6.05 

4.85 

4.52 

4.46 

4.46 

4.44 

4.39 

4.36 

4.29 

4.14 

3.95 

0.8 

12.24 

6.90 

5.36 

4.93 

4.85 

4.85 

4.83 

4.77 

4.74 

4.66 

4.50 

4.29 

0.9 

14.26 

7.75 

5.85 

5.30 

5.19 

5.19 

5.18 

5.11 

5.08 

4.99 

4.82 

4.59 

1.0 

16.36 

8.61 

6.32 

5.64 

5.50 

5.51 

5.49 

5.42 

5.38 

5.29 

5.11 

4.86 

1.1 

18.53 

9.47 

6.77 

5.95 

5.77 

5.90 

5.87 

5.70 

5.66 

5.56 

5.37 

5.10 

1.2 

20.77 

10.34 

7.21 

6.23 

6.01 

6.24 

6.26 

5.95 

5.90 

5.80 

5.59 

5.32 

1.3 

23.07 

11.22 

7.63 

6.50 

6.23 

6.54 

6.61 

6.24 

6.12 

6.01 

5.79 

5.50 

1.4 

25.43 

12.10 

8.04 

6.74 

6.42 

6.81 

6.93 

6.58 

6.36 

6.20 

5.97 

5.67 

1.5 

27.84 

12.99 

8.45 

6.97 

6.59 

7.05 

7.21 

6.89 

6.67 

6.37 

6.13 

5.81 

1.6 

30.30 

13.89 

8.84 

7.18 

6.75 

7.25 

7.46 

7.16 

6.95 

6.52 

6.27 

5.94 

1.7 

32.81 

14.79 

9.23 

7.38 

6.88 

7.44 

7.69 

7.41 

7.20 

6.66 

6.39 

6.04 

1.8 

35.36 

15.71 

9.61 

7.56 

7.00 

7.59 

7.89 

7.64 

7.43 

6.85 

6.49 

6.14 

1.9 

37.96 

16.63 

9.98 

7.74 

7.11 

7.73 

8.07 

7.84 

7.63 

7.06 

6.58 

6.21 

2.0 

40.59 

17.55 

10.36 

7.90 

7.20 

7.85! 

8.23 

8.02 

7.82 

7.24 

6.66 

6.28 
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15%  Diffuser  Loss  Data: 


First  Solution: 


Table  20:  Table  of  Entropy  Change  15%  Diffuser 
Loss  (1st  Solution) 


Values  of 

ASinr 

Pr  =  0.9003 

Pr  =  0.10 

Pr  =  0.55 

Pr  =  0.4 

rtii  R 
//=0.0 
/y=  1.0 

/i  =  2.0 

999E-09 

0.980067 

1.423613 

3.33E-'* 

0.890117 

1.312721 

1.44E'‘" 

0.895577 

1.282924 

-3.9E‘'^ 

0.986324 

1.394148 

Table  21:  Thrust  Augmentation  with  15% 
Diffuser  Loss  (1st  Solution) 


Thrust  when  expanded  to  ambient  Pressure  - 
%  Thrust  increase  over  baseline,  -107  kN,  Diffuser  loss=15% 


M//Pr= 

0.90 

0.85 

0.80 

0.75 

0.70 

0.65 

0.60 

0.55 

0.5328 

0.50 

0.45 

0.40 

0.39 

0.49 

0.57 

0.65 

0.73 

0.81 

0.89 

0.97 

1.0000 

1.05 

1.13 

1.20 

0.0 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.1 

0.12 

0.25 

0.26 

0.18 

0.03 

-0.21 

-0.53 

-0.96 

-1.14 

-1.52 

-222 

-3.10 

0.2 

0.45 

0.76 

0.89 

0.91 

0.83 

0.66 

0.39 

0.00 

-0.16 

-0.52 

-1.20 

-2.09 

0.3 

0.77 

1.22 

1.44 

1.53 

1.51 

1.40 

1.17 

0.83 

0.67 

0.33 

-0.34 

-1.23 

0.4 

1.06 

1.62 

1.91 

2.05 

2.09 

2.02 

1.84 

1.52 

1.38 

1.05 

0.39 

-0.50 

0.5 

1.33 

1.97 

2.31 

2.50 

2.57 

2.54 

2.40 

2.12 

1.98 

1.67 

1.02 

0.13 

0.6 

1.57 

2.28 

2.66 

2.88 

2.99 

2.99 

2.88 

2.62 

2.50 

2.20 

1.56 

0.66 

0.7 

1.78 

2.54 

2.96 

3.21 

3.34 

3.37 

3.28 

3.05 

2.93 

2.65 

2.02 

1.12 

0.8 

1.97 

2.77 

3.22 

3.49 

3.64 

3.68 

3.62 

3.42 

3.31 

3.03 

2.41 

1.51 

0.9 

2.12 

2.96 

3.44 

3.72 

3.89 

3.95 

3.91 

3.72 

3.62 

3.36 

2.75 

1.84 

1.0 

2.25 

3.12 

3.62 

3.92 

4.09 

4.17 

4.14 

3.98 

3.89 

3.63 

3.03 

2.12 

1.1 

2.36 

3.26 

3.78 

4.09 

4.27 

4.35 

4.34 

4.20 

4.11 

3.86 

3.27 

2.36 

1.2 

2.45 

3.37 

3.90 

4.22 

4.41 

4.49 

4.49 

4.37 

4.29 

4.05 

3.47 

2.55 

1.3 

2.51 

3.45 

4.00 

4.33 

4.52 

4.61 

4.64 

4.51 

4.43 

4.21 

3.64 

2.72 

1.4 

2.56 

3.52 

4.08 

4.41 

4.60 

4.73 

4.84 

4.62 

4.55 

4.34 

3.77 

2.84 

1.5 

2.59 

3.56 

4.13 

4.47 

4.67 

4.86 

5.02 

4.70 

4.64 

4.43 

3.88 

2.94 

1.6 

2.60 

3.59 

4.17 

4.52 

4.71 

4.96 

5.17 

4.86 

4.70 

4.50 

3.96 

3.02 

1.7 

2.60 

3.60 

4.19 

4.54 

4.74 

5.04 

5.29 

5.00 

4.78 

4.55 

4.01 

3.07 

1.8 

2.59 

3.60 

4.19 

4.55 

4.75 

5.09 

5.39 

5.12 

4.91 

4.58 

4.05 

3.10 

1.9 

2.56 

3.58 

4.18 

4.54 

4.74 

5.12 

5.46 

5.22 

5.01 

4.59 

4.06 

3.11 

2.0 

2.52 

3.55 

4.15 

4.52 

4.72 

5.13 

5.52 

5.30 

5.09 

4.58 

4.06 

3.11 
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Second  Solution: 


Table  22:  Table  of  Entropy  Change  15%  Diffuser 
Loss  (2nd  Solution) 


Values  of 

ASinl 

Fr  =  0.9003 

Fr  =  0.70 

Fr  =  0.55 

Fr  =  0.4 

• 

mA  R 

/j=0.0 

THE*'" 

3.33E-'7 

1.44E''- 

-3.9E‘'^ 

M=  1.0 

0.30201 

0.883587 

0.887248 

0.912787 

//=2.0 

-0.67188 

1.305672 

1.282924 

1.351426 

Table  23:  Thrust  Augmentation  with  1 5% 
Diffuser  Loss  (2nd  Solution) 


Thrust  when  expanded  to  ambient  Pressure  • 

7o  Thrust  increase  over  baseline,  MOT  kN,  Diffuser  loss=15®/o 


VJIPr^ 

0.90 

0.85 

0.80 

0.75 

0.70 

0.65 

0.60 

0.55 

0.5328 

0.50 

0.45 

0.40 

0.39 

0.49 

0.57 

0.65 

0.73 

0.81 

0.89 

0.97 

1.0000 

1.05 

1.13 

1,20 

0.0 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.1 

0.83 

0.75 

0.74 

0.75 

0.76 

0.76 

0.75 

0.74 

0.73 

0.72 

0.69 

0.65 

0.2 

1.89 

1.49 

1.41 

1.40 

1.41 

1.41 

1.40 

1.37 

1.36 

1.34 

1.28 

1.21 

0.3 

3.14 

2.23 

2.01 

1.97 

1.97 

1.97 

1.96 

1.93 

1.91 

1.87 

1.80 

1.70 

0.4 

4.55 

2.97 

2.56 

2.46 

2.45 

2.45 

2.44 

2.40 

2.38 

2.34 

2.25 

2.12 

0.5 

6.09 

3.70 

3.06 

2.90 

2.87 

2.87 

2.86 

2.81 

2.79 

2.74 

2.63 

2.48 

0.6 

7.75 

4.44 

3.53 

3.28 

3.23 

3.23 

3.22 

3.17 

3.15 

3.09 

2.96 

0.7 

9.51 

5.19 

3.96 

3.61 

3.54 

3.54 

3.52 

3.48 

3.45 

3.38 

3.25 

3.06 

0.8 

11.38 

5.94 

4.36 

3.91 

3.81 

3.80 

3.79 

3.74 

3.71 

3.64 

3.49 

3.29 

0.9 

13.33 

6.69 

4.75 

4.16 

4.03 

4.03 

4.02 

3.96 

3.Q3 

3.86 

3.70 

3.48 

1.0 

15.36. 

7.45 

5.11 

4.39 

4.23 

4.22 

4.21 

4.15 

4.12 

4.04 

3.87 

3.63 

1.1 

17.47 

8.22 

5.46 

4.59 

4.39 

4.37 

4.37 

4.31 

4.28 

4.19 

4.01 

3.76 

1.2 

19.65 

9.00 

5.79 

4.77 

4.52 

4.50 

4.50 

4.44 

4.41 

4.32 

4.12 

1.3 

21.89 

9.79 

6.11 

4.93 

4.63 

4.61 

4.64 

4.55 

4.51 

4.42 

4.21 

3.94 

1.4 

24.19 

10.58 

6.42 

5.07 

4.71 

4.73 

4.84 

4.63 

4.59 

4.49 

4.28 

3.99 

1.5 

26.55 

11.39 

6.73 

5.19 

4.77 

4.86 

5.02 

4.70 

4.65 

4.55 

4.33 

4.02 

1.6 

28.96 

12.20 

7.02 

5.30 

4.82 

4.96 

5.17 

4.86 

4.70 

4.59 

4.36 

4.04 

1.7 

31.42 

13.02 

7.31 

5.39 

4.85 

5.04 

5.29 

5.00 

4.78 

4.61 

4.37 

4.04 

1.8 

33.93 

13.85 

7.60 

5.48 

4.87 

5.09 

5.39 

5.12 

4.91 

4.61 

4.36 

4.03 

1.9 

36.48 

14.68 

7.88 

5.55 

4.87 

5.12 

5.46 

5.22 

5.01 

4.61 

4.34 

4.00 

2.0 

39.07 

15.53 

8.15 

5.61 

4.85 

5.13 

5.521 

5.30 

5.09 

4.58 

4.31 

3.95 
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20%  Diffuser  Loss  Data: 


First  Solution: 


Table  24:  Table  of  Entropy  Change  20%  Diffuser 
Loss  (1st  Solution) 


Values  of 

ASioi 

Pr  =  0.9003 

Pr  =  0.70 

Pr  =  0.55 

Pr  =  0.4 

• 

nii  R 

//=0.0 

1.12E'°* 

-1.2E''" 

-1.7E''^ 

-6.1E-‘* 

//=  1.0 

1.004547 

0.903166 

0.909141 

1.004222 

II 

b 

1.451494 

1.325713 

1.311922 

1.409884 

Table  25:  Thrust  Augmentation  with  20% 
Diffuser  Loss  ( 1  st  Solution) 


Thrust  when  expanded  to  ambient  Pressure  • 

%  Thrust  increase  over  baseline,  ~107  kN,  Diffuser  ioss=20% 


vJIPr^ 

0.90 

0.85 

0.80 

0.75 

0.70 

0.65 

0.60 

0.55 

0.5328 

0.50 

0.45 

0.40 

0.39 

0.49 

0.57 

0.65 

0.73 

0.81 

0.89 

0.97 

1.0000 

1.05 

1.13 

1.20 

0.0 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.1 

-0.24 

-0.09 

-0.09 

-0.19 

-0.36 

-0.63 

-0.99 

-1.45 

-1.64 

-2.04 

-2.78 

-3.71 

0.2 

-0.11 

0.24 

0.38 

0.40 

0.31 

0.12 

-0.18 

-0.60 

-0.78 

-1,16 

-1.88 

-2.82 

0.3 

0.03 

0.53 

0.78 

0.88 

0.86 

0.73 

0.48 

0.10 

-0.06 

-0.43 

-1.14 

-2.07 

0.4 

0.16 

0.78 

1.11 

1.27 

1.31 

1.23 

1.02 

0.68 

0.53 

0.18 

-0.52 

-1.45 

0.5 

0.27 

0.98 

1.38 

1.59 

1.67 

1.63 

1.47 

1.16 

1.02 

0.68 

0.00 

-0.94 

0.6 

0.35 

1.15 

1.59 

1.84 

1.96 

1.96 

1.83 

1.55 

1.42 

1.09 

0.42 

-0.52 

0.7 

0.42 

1.27 

1.76 

2.04 

2.19 

2.22 

2.12 

1.87 

1.74 

1.43 

0.77 

-0.18 

0.8 

0.46 

1.37 

1.88 

2.20 

2.37 

2.42 

2.35 

2.12 

2.00 

1.70 

1.05 

0.10 

0.9 

0.48 

1.43 

1.97 

2.31 

2.50 

2.57 

2.52 

2.32 

2.20 

1.92 

1.27 

0.32 

1.0 

0.48 

1.46 

2.03 

2.38 

2.59 

2.68 

2.65 

2.46 

2.36 

2.08 

1.45 

0.49 

1.1 

0.46 

1.47 

2.06 

2.43 

2.64 

2.75 

2.73 

2.57 

2.47 

2.20 

1.57 

0.61 

1.2 

0.41 

1.45 

2.06 

2.44 

2.67 

2.78 

2.78 

2.63 

2.54 

2.28 

1.66 

0.69 

1.3 

0.35 

1.41 

2.04 

2.43 

2.66 

2.78 

2.79 

2.67 

2.58 

2.33 

1.72 

0.74 

1.4 

0.28 

1.36 

1.99 

2.39 

2.63 

2.76 

2.78 

2.67 

2.58 

2.35 

1.74 

0.76 

1.5 

0.18 

1.28 

1.93 

2.33 

2.58 

2.71 

2.74 

2.64 

2.56 

2.33 

1.73 

0.75 

1.6 

0.07 

1.19 

1.85 

2.26 

2.51 

2.64 

2.67 

2.59 

2.52 

2.30 

1.70 

0.71 

1.7 

-0.05 

1.08 

1.75 

2.17 

2.42 

2.55 

2.67 

2.52 

2,45 

2.24 

1.65 

0.65 

1.8 

-0.18 

0.95 

1.63 

2.06 

2.31 

2.44 

2.65 

2.42 

2.36 

2.16 

1.58 

0.57 

1.9 

-0.33 

0.82 

1.50 

1.93 

2.19 

2.32 

2.61 

2.36 

2.25 

2.06 

1.48 

0.48 

2.0 

-0.49 

0.67 

1.36 

■  1.79 

2.05 

2:i8 

2.56 

2.33 

2.13 

1.94 

1.37 

0.36 
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Second  Solution: 


Table  26:  Table  of  Entropy  Change  20%  Diffuser 
Loss  (2nd  Solution) 


Values  of 

ASfoi 

Pr  =  0.9003 

P.  =  0.70 

Pr  =  0.55 

Pr  =  0.4 

m^  R 

//=0.0 

2.78E‘‘^ 

-1.2E'‘-" 

-1.7E-'-' 

-6.1E-'^ 

H=  1.0 

0.29853 

0.890359 

0.893876 

0.917363 

//=2.0 

-0.69009 

1.312635 

1.311922 

1.356806 

Table  27:  Thrust  Augmentation  with  20% 
Diffuser  Loss  (2nd  Solution) 


Thrust  when  expanded  to  ambient  Pressure  - 
%  Thrust  increase  over  baseline,  -107  kN,  Diffuser  loss=:20% 


M//Pr= 

0.90 

0.85 

0.80 

0.75 

0.70 

0.65 

0.60 

0.55 

0.5328 

0.50 

0.45 

0.40 

0.39 

0.49 

0.57 

0.65 

0.73 

0.81 

0.89 

0.97 

1.0000 

1.05 

1.13 

1.20 

0.0 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.1 

0.66 

0.58 

0.58 

0.59 

0.60 

0.60 

0.60 

0.59 

0.58 

0.57 

0.54 

0.50 

0.2 

1.58 

1.18 

1.10 

1.09 

1.10 

1.10 

1.10 

1.08 

1.07 

1.04 

0.99 

0.92 

0.3 

2.71 

1.78 

1.56 

1.51 

1.51 

1.52 

1.51 

1.48 

1.47 

1.44 

1.37 

1.27 

0.4 

4.01 

2.39 

1.97 

1.87 

1.86 

1.86 

1.85 

1.82 

1.80 

1.76 

1.68 

1.55 

0.5 

5.45 

3.01 

2.34 

2.17 

2.14 

2.14 

2.13 

2.09 

2.08 

2.03 

1.93 

1.78 

0.6 

7.01 

3.63 

2.68 

2.42 

2.37 

2.36 

2.35 

2.32 

2.29 

2.24 

2.13 

1.97 

0.7 

8.70 

4.26 

2.99 

2.63 

2.55 

2.54 

2.53 

2.49 

2.47 

2.41 

2.28 

2.10 

0.8 

10.48 

4.90 

3.28 

2.80 

2.68 

2.67 

2.66 

2.62 

2.60 

2.53 

2.40 

2.20 

0.9 

12.36 

5.55 

3.55 

2.93 

2.78 

2.77 

2.76 

2.72 

2.69 

2.62 

2.47 

2.26 

1.0 

14.32 

6.21 

3.80 

3.04 

2.85 

2.83 

2.82 

2.78 

2.75 

2.68 

2.52 

2.30 

1.1 

16.37 

6.89 

4.03 

3.13 

2.89 

2.86 

2.86 

2.81 

2.78 

2.71 

2.54 

2.30 

1.2 

18.48 

7.57 

4.25 

3.19 

2.90 

2.87 

2.86 

2.82 

2.79 

2.71 

2.53 

2.28 

1.3 

20.66 

8.26 

4.47 

3.23 

2.89 

2.85 

2.85 

2.80 

2.77 

2.69 

2.50 

2.23 

1.4 

22.91 

8.96 

4.67 

3.25 

2.86 

2.80 

2.81 

2.76 

2.73 

2.64 

2.45 

2.17 

1.5 

25.21 

9.67 

4.87 

3.26 

2.80 

2.74 

2.75 

2.70 

2.67 

2.58 

2.37 

2.08 

1.6 

27.57 

10.39 

5.06 

3.26 

2.74 

2.66 

2.67 

2.62 

2.59 

2.49 

2.28 

1.98 

1.7 

29.98 

11.11 

5.24 

3.24 

2.65 

2.56 

2.67 

2.53 

2.49 

2.39 

2.18 

1.87 

1.8 

32.44 

11.85 

5.42 

3.22 

2.55 

2.45 

2.65 

2.42 

2.38 

2.28 

2.05 

1.73 

1.9 

34.95 

12.60 

5.60 

3.18 

2.44 

2.32 

2.61 

2.36 

2.26 

2.15 

1.92 

1.59 

2.0 

37.49 

13.36 

5.77 

3.13 

2.32 

2.18 

2.56 

2.33 

2.13 

2.01 

1.77 

1.43 
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25%  Diffuser  Loss  Data: 


First  Solution: 


Table  28:  Table  of  Entropy  Change  25%  Diffuser 
Loss  (1st  Solution) 


Values  of 

ASioi 

=  0.9003 

P,  =  0.70 

Pr  =  0.55 

II 

o 

m 

mi  R 

2.78E'”^ 

3.33E''^ 

5.55E''^ 

-2.3E'"' 

1.0 

1.032247 

0.919235 

0.925613 

1.024761 

/i=2.0 

1.482709 

1.341759 

1.332471 

1.428351 

Table  29:  Thrust  Augmentation  with  25% 
Diffuser  Loss  (1st  Solution) 


Thrust  when  expanded  to  ambient  Pressure  - 
%  Thrust  increase  over  baseline,  --107  kN,  Diffuser  loss=25% 


M//Pr= 

0.90 

0.85 

0.80 

0.75 

0.70 

0.65 

0.60 

0.55 

0.5328 

0.50 

0.45 

0.40 

0.39 

0.49 

0.57 

0.65 

0.73 

0.81 

0.89 

0.97 

1.0000 

1.05 

1.13 

1.20 

0.0 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.1 

-0.66 

-0.50 

-0.51 

-0.62 

-0.82 

-1.11 

-1.50 

-2.00 

-2.20 

-2.63 

-3.41 

-4.39 

0.2 

-0.75 

-0.35 

-0.20 

-0.19 

-0.29 

-0.51 

-0.83 

-1.28 

-1.47 

-1.88 

-2.64 

-3.62 

0.3 

-0.81 

-0.24 

0.04 

0.14 

0.11 

-0.03 

-0.30 

-0.71 

-0.89 

-1.28 

-2-03 

-3.00 

0.4 

-0.86 

-0.16 

0.21 

0.39 

0.42 

0.33 

0.11 

-0.26 

-0.42 

-0.80 

-1.53 

-2.51 

0.5 

-0.92 

-0.12 

0.32 

0.56 

0.65 

0.61 

0.43 

0.09 

-0.06 

-0.42 

-1.14 

-2.13 

0.6 

-1.00 

-0.11 

0.39 

0.68 

0.81 

0.81 

0.66 

0.36 

0.21 

-0.13 

-0.84 

-1.83 

0.7 

-1.09 

-0.13 

0.42 

0.74 

0.91 

0.94 

0.83 

0.55 

0.42 

0.08 

-0.62 

-1.61 

0.8 

-1.20 

-0.19 

0.40 

0.76 

0.95 

1.01 

0.93 

0.68 

0.55 

0.23 

-0.46 

-1.46 

0.9 

-1.33 

-0.27 

0.35 

0.73 

0.95 

1.04 

0.98 

0.76 

0.63 

0.33 

-0-36 

-1.36 

1.0 

-1.48 

-0.38 

0.27 

0.68 

0.92 

1.02 

0.98 

0.78 

0.67 

0.37 

-0.31 

-1.32 

1.1 

-1.64 

-0.51 

0.16 

0.59 

0.84 

0.96 

0.95 

0.76 

0.65 

0.37 

-0.30 

-1.31 

1.2 

-1.82 

-0.66 

0.03 

0.47 

0.74 

0.87 

0.87 

0.71 

0.61 

0.33 

-0.33 

-1.35 

1.3 

-2.02 

-0.84 

-0.12 

0.33 

0.61 

0.75 

0.77 

0.62 

0.52 

0.25 

-0.40 

-1 .43 

1.4 

-2.23 

-1.03 

-0.30 

0.16 

0.45 

0.61 

0.63 

0.50 

0.41 

0.15 

-0.50 

-1.53 

1.5 

-2.46 

-1.24 

-0.50 

-0.02 

0.27 

0.43 

0.47 

0.36 

0.27 

0.02 

-0.63 

-1.66 

1.6 

-2.70 

-1.46 

-0.71 

-0.23 

0.07 

0.24 

0.29 

0.19 

0.11 

-0.14 

-0.78 

-1.82 

1.7 

-2.96 

-1.70 

-0.94 

-0.45 

-0.14 

0.03 

0.09 

0.00 

-0.08 

-0.32 

-0.95 

-2.00 

1.8 

-3.22 

-1.95 

-1.18 

-0.68 

-0.37 

-0.20 

-0.14 

-0.21 

-0.29 

-0.52 

-1.14 

-2.20 

1.9 

-3.50 

-2.22 

-1.44 

-0.94 

-0.62 

-0.44 

-0.38 

-0.44 

-0.51 

-0.74 

-1.36 

-2.42 

2.0 

-3.78 

-2.49 

-1.70 

-1.20 

-0.88 

-0.70 

-0.63 

-0.69 

-0.76 

-0.97 

-1.59 

-2.66 

Second  Solution: 


Table  30:  Table  of  Entropy  Change  25%  Diffuser 
Loss  (2nd  Solution) 


Values  of 

A.Slnt 

Pr  =  0.9003 

Pr  =  0.70 

Pr  =  0.55 

II 

o 

m-i  R 
//  =  0.0 

1.0 
/r  =  2.0 

3.33E*'^ 

0.29339 

-0.71263 

3.33E*'^ 

0.898159 

1.320856 

5.55E''® 

0.901529 

1.331072 

-2.3E‘‘-" 

0.92301 

1.363542 

Table  31 :  Thrust  Augmentation  with  25% 
Diffuser  Loss  (2nd  Solution) 


Thrust  when  expanded  to  ambient  Pressure  - 


vJfPr^ 

0.90 

%  Thrust  increase  over  baseline,  -107  kN,  Diffuser  loss: 
0.85  0.80  0.75  0.70  0.65  0.60 

=25% 

0.55 

0.5328 

0.50 

0.45 

0.40 

0.39 

0.49 

0.57 

0.65 

0.73 

0.81 

0.89 

0.97 

1.0000 

1.05 

1.13 

1.20 

0.0 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.1 

0.48 

0.40 

0-40 

0.42 

0.43 

0.43 

0.43 

0.42 

0-41 

0.40 

0.38 

0.34 

0.2 

1.25 

0.84 

0.76 

0.75 

0.76 

0.77 

0.77 

0.75 

0.74 

0.72 

0.67 

0.60 

0.3 

2.25 

1.29 

1.06 

1.02 

1.02 

1.03 

1.02 

1.00 

0.99 

0.96 

0.89 

0.80 

0.4 

3.43 

1.77 

1.33 

1.22 

1.21 

1.22 

1.21 

1.18 

1.17 

1.13 

1.05 

0.94 

0.5 

4.77 

2.25 

1.56 

1.37 

1.34 

1.34 

1.34 

1.31 

1.29 

1.25 

1.16 

1.02 

0.6 

6.24 

2.75 

1.77 

1.48 

1.42 

1.42 

1.42 

1.39 

1,37 

1.32 

1.21 

1.06 

0.7 

7.83 

3.27 

1.94 

1.55 

1.46 

1.45 

1.45 

1.42 

1.40 

1.34 

1.23 

1.05 

0.8 

9.54 

3.79 

2.10 

1.59 

1.46 

1.44 

1.44 

1.41 

1.38 

1.33 

1.20 

1.01 

0.9 

11.34 

4.33 

2.24 

1.59 

1.42 

1.40 

1.40 

1.36 

1.34 

1.28 

1.14 

0.94 

1.0 

13.23 

4.88 

2.37 

1.57 

1.36 

1.32 

1.32 

1.29 

1.26 

1.20 

1.05 

0.84 

1.1 

15.20 

5.44 

2.48 

1.53 

1.26 

1.22 

1.22 

1.18 

1.16 

1.09 

0.93 

0.71 

1.2 

17.26 

6.02 

2.59 

1.46 

1.14 

1.09 

1.09 

1.05 

1.03 

0.96 

0.79 

0.55 

1.3 

19.38 

6.60 

2.68 

1.38 

1.00 

0.94 

0.93 

0.90 

0.87 

0.80 

0.63 

0.38 

1.4 

21.56 

7.20 

2.77 

1.28 

0.84 

0.76 

0.76 

0.72 

0.70 

0.62 

0.44 

0.18 

1.5 

23.81 

7.81 

2.85 

1.16 

0.66 

0.57 

0.57 

0.53 

0.50 

0.43 

0.24 

-0.03 

1.6 

26.12 

8.43 

2.92 

1.04 

0.46 

0.36 

0.36 

0.32 

0.29 

0.21 

0.02 

-0.26 

1.7 

28.48 

9.06 

2.99 

0.90 

0.25 

0.13 

0.13 

0.10 

0.07 

-0.02 

-0.21 

-0.51 

1.8 

30.89 

9.70 

3.06 

0.75 

0.03 

-0.11 

-0.11 

-0.14 

-0.17 

-0.26 

-0.46 

-0.77 

1.9 

33.34 

10.35 

3.12 

0.59 

-0.21 

-0.37 

-0.36 

-0.40 

-0.43 

-0.52 

-0.73 

-1.04 

2.0 

35.84 

11.02 

3.18 

0.43 

-0.45 

-0.64 

-0.63 

-0.66 

-0.69 

-0.79 

-1.00 

-1.32 

